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ABSTRACT 



14 The LAT instrument on the Fermi mission will reveal the rich spectral and 
temporal gamma-ray burst phenomena in the >100 MeV band. The synergy 
with Fermi's GBM detectors will link these observations to those in the well- 
explored 10-1000 keV range; the addition of the >100 MeV band observations 
will resolve theoretical uncertainties about burst emission in both the prompt 
and afterglow phases. Trigger algorithms will be applied to the LAT data both 
onboard the spacecraft and on the ground. The sensitivity of these triggers will 
differ because of the available computing resources onboard and on the ground. 
Here we present the LAT's burst detection methodologies and the instrument's 
GRB capabihties. 

15 Subject headings: gamma rays: bursts 



16 1. Introduction 

17 The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope (formerly 

18 GLAST — Gamma-ray Large Area Space Telescope) will turn the study of the 20 MeV to more 

19 than 300 GeV spectral and temporal behavior of gamma-ray bursts (GRBs) from specula- 

20 tion based on a few suggestive observations to a decisive diagnostic of the emission processes. 

21 The burst observations of the Energetic Gamma-Ray Experiment Telescope (EGRET) on 
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22 the Compton Gamma-Ray Observatory (CGRO) suggested three types of high energy emis- 

23 sion: an extrapolation of the 10-1000 keV spectral component to the >100 MeV band; an 

24 additional spectral component during the <1 MeV 'prompt' emission; and high energy emis- 

25 sion that lingers long after the prompt emission has faded away. The LAT's observations, 

26 in conjunction with the Gamma-ray Burst Monitor (GBM — 8 keV to 30 MeV), will provide 

27 unprecedented spectral-temporal coverage for a large number of bursts. The spectra from 
2B these two instruments will cover seven and a half energy decades (<10 keV to >300 GeV; 
29 see Fig. [H which shows different theoretically- predicted spectra). Thus the LAT will explore 

the rich phenomena suggested by the EGRET observations, probing the physical processes 



30 



31 in the extreme radiating regions. 

32 In this paper we provide the scientific community interested in GRBs with an overview of 

33 the LAT's operations and capabilities in this research area. Our development of detection and 

34 analysis tools has been guided by the previous observations and the theoretical expectations 

35 for e mission in the >100 M eV band (§ [2]). The LAT is described in depth in an instrument 



36 paper (lAtwood et al.ll2009l ). and therefore here we only provide a brief summary of the Fermi 

37 mission and the LAT, focusing on issues relevant to burst detection and analysis (§ |3]). 
3B Simulations are the basis of our analysis of the mission's burst sensitivity, and are largely 

39 based on CGRO observations (§ H]). We use our simulation methodology to estimate the 

40 ultimate burst sensitivity and the resulting burst flux distribution (§|5]). Both the LAT and 

41 the GBM will apply burst detection algorithms onboard and on the ground, and the efficiency 

42 of these methods will determine which bursts the LAT will detect, and with what latency 

43 (§ ED- Once a burst has been detected, spectral and temporal analysis of LAT (and GBM) 

44 data will be possible (§ [7]). The burst observations by ground-based telescopes and other 

45 space missions, particularly Swift, will complement the Fermi observations (§ |8]). While 

46 basic methods are in place for detecting and analyzing burst data, in-fiight experience will 

47 guide future work (§[n])- 
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Fig. 1. — Simulated gamma-ray burst spectra, showing the broad energy range covered by 
Fermi: (from left to right) the GBM Nal (blue band: 8-2000 keV), the GBM BGO (green: 
150 keV-30 MeV) and the LAT (red curve: 20 MeV to >300 GeV) detectors. The dashed 
curves are simple extrapolations of the typical GRB 10-1000 keV spectra into the GeV band, 
while the solid curves add an exponential cutoff that might result from absorption internal or 
external to the burst. The two different high energy photon indices (3=-2.25 (black curves) 
and P=-2.5 (grey curves) demonstrate the dependence of the expected LAT flux on this 
photon index. There may be additional high energy components that are not known yet and 
are not shown in the figure. 
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2. Burst Physics Above 100 MeV 



2.1. Previous Observations 
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The detectors of the Compton Gamma-Ray Observatory ( CGRO) provided time- resolved 
spectra for a statistically well-defined burst population. These observations are the foun- 
dation of our expectations for FermVs discoveries, which have guided the development of 
analysis tools before launch. 

The Burst And Transient Source Experiment (BATSE) on CGRO observed a large 
sample of bursts in the ~25-2000 keV band with well-understood population statistics 
(jPaciesas et al.lll999l ). Spectroscopy by the BATSE detectors found that the emi ssion in this 



energy band could be described by the empirical four parameter "Band" function ( iBand et al. 
1993h 



E°exp[-E(2 + a)/Ep] 



NBe.nd{E\No,Ep,a,p) = No 



where a and P are the low and high energy photon indices, respectively, and Ep is the 'peak 
energy' which corresponds to the maximum of E'^N( E) oc uf,, for the low energy compo- 



nent. Typically a ^ —0.5 to —1 and f5 is less than —2 (iBand et al.lll993l : iPreece et al.l 12000 



Kaneko et al.ll2006l ): the total energy would be infinite if /9 > —2 unless the spectrum has 
a high energy cutoff. The observations of 37 bursts by the Compton Telescope (COMP- 
TEL) on C GRO (0.75-30 MeV) are consistent with the BATSE observations of this spectral 
component ( Hoover et al.| 2005 ). Because of the relatively poor spectral resolution of the 
BATSE detectors (lBriggslll999l ). this functional form usually is a good description of spec- 
tra accumulated over bot h short tirne per iods and entire bursts, even though bursts show 
strong spectral evolution (IFord et al.lll995l ). It is this 10-1000 keV 'prompt' component that 
is well-characterized and therefore provides a basis for quantitative predictions. A detailed 
duration-integrated spectral analysi s (in 30 keV- 2 00 Me V) of the prompt emission for 15 
bright BATSE GRB performed by iKaneko et al.l (l2008l ) confirmed that only in few case 
there's a significant high-energy excess with respect to low energy spectral extrapolations. 

The burst observations by the Energetic Gamma-Ray Experiment Telescope (EGRET) 
on CGRO (20 MeV to 30 GeV) provide the best prediction of the LAT observations. EGRET 
observed different types of high energy burst phenomena. Four bursts had simultaneous 
emission in both the EGRET and BATSE energy bands, su ggesting that the spectrum 
observed by BATSE extrapolates to the EGRET energy band (jDingud 120031 ). However, the 
correlation with the prompt phase pulses was hampered by the severe EGRET spark chamber 
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dead time (~100 ms/event) that was comparable or longer than the pulse timescales. The 
EGRET observations of these bursts suggest that the ~1 GeV emission often lasts longer 
than the lower energy emission, and thus results in part from a differ ent physical origin. A 
similar behaviour is present also in GRB 080514B detected by AGILE (jGiuliani et al.ll2008l ). 



Whether high energy emission is present in both long and short bursts is unknown. The 
four bursts with high energy emission detected by EGRET were all long bursts, alth ough 
GRB 930131 is an interesting case. It was detected by BATSE (IKouveliotou et al.lll994l ) with 
duration of T9o=14 s0 and found to have high- energy (>30 MeV) p hotons accompanying the 
prompt phase and possibly extending beyond (jSommer et al.lll994l ). The BATSE lightcurve 
is dominated by a hard initial emission lasting 1 sec and followed by a smooth extended 
emission. This burst may, therefore, have been one of those long bursts possibly associated 
with a merger and not a collapsar ori gin, common ly understood as the most probable origin 
for short and long burst respecti velvdZhane] 120071) . Sever al events have now been identified 
that could fit into this category (INorris fc Bonnelll l2006l ) and their origin is still uncertain. 
LAT will make an important contribution in determining the nature of the high energy 
emission from similar events and a larger sample of bursts with detected high energy emission 
will determine whether the absence of high energy emission differentiates short from long 
bursts. 

A high energy temporally res olved spectral compon ent in addition to the Band function 
is clearly present in GRB 941017 (iGonzalez et al.ll2003l ): this component is harder than the 
low energy prompt component, and continues after the low energy component fades into the 
background. The time integrat ed spectra of both G RB 941017 and GRB 980923 show this 
additional spectral component (IKaneko et al.ll2008l ). 



Finally, the >1 GeV emission lingered for 90 minutes after the prompt low energy 
emission for GRB 940217, including an 18 GeV photon 1.5 hours after the burst trigger 
f lHurley et al.lll994l ). Whether this emission is physically associated with the lower energy 
afterglows is unknown. 

These three empirical types of high energy emission — an extrapolation of the low energy 
spectra; an additional spectral component during the low energy prompt emission; and an 
afterglow — guide us in evaluating Fermz's burst observation capabilities. 

Because the prompt low energy component was characterized quantitatively by the 
BATSE observations while the EGRET observations merely demonstrated that different 
components were present, our simulations are based primarily on extrapolations of the 



^Tgo is the time over which 90% of the emission occurs in a specific energy band. 
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112 prompt low energy component from the BATSE band to the >100 MeV band. We rec- 

113 ognize that the LAT will probably detect additional spectral and temporal components, or 

114 spectral cutoffs, that are not treated in this extrapolation. 

115 During the first few months of the Ferm i mission, L AT detected already emission from 

116 three GRBs: 080825C jBouvier et all l2008h . 080916C JXaiima et all hoOS^ and 081024B 

117 (jOmodeil |2008| ) . The rich phenomenology of high energy emission is confirmed in these 
iiB three events, where spectral measurements over various order of magnitude were possible 

119 together with the detection of extended emission and spectral lags. In particular, the GRB 

120 080916C was bright e nough to afford u nprecedented broad-band spectral coverage in four 

121 distinct time intervals (lAbdo et al.ll2009l ). thereby offering new insights into the character of 

122 energetic bursts. 



123 2.2. Theoretical Expectations 

124 In the current standard scenario, the burst emission arises in a highly relativistic, un- 

125 steady outflow. Several different progenitor types could create this outflow, but the initial 

126 high optical depth within the outflow obscures the progenitor type. As this outflow gradu- 

127 ally becomes optically thin, dissipation processes within the outflow, as well as interactions 
12B with the surrounding medium, cause particles to be accelerated to high energies and loose 

129 some of their energy into radiation. Magnetic fields at the emission site can be strong and 

130 may be caused by a frozen-in component carried out by the outfiow from the progenitor, or 

131 may be built up by turbulence or collisionless shocks. The emitted spectral distribution then 

132 depends on the details of the radiation mechanism, particle acceleration, and the dynamics 

133 of the explosion itself. 

134 'Internal shocks' result when a faster region catches up with a slower region within 

135 the outfiow. 'External shocks' occur at the interface between the outfiow and the ambient 

136 medium, and include a long-lived forward shock that is driven into the external medium 

137 and a short-lived reverse shock that decelerates the outfiow. Thus the simple model of a 
13B one-dimensional relativistic outfiow leads to a multiplicity of shock fronts, and many possible 

139 interacting emission regions. 

140 As a result of the limited energy ranges of past and current experiments, most theories 

141 have not been clearly and unambiguously tested. Fermi's GBM and LAT will provide 

142 an energy range broad enough to distinguish between different origins of the emission; in 

143 particular the unprecedented high-energy spectral coverage will constrain the total energy 

144 budget and radiative efficiency, as potentially most of the energy may be radiated in the LAT 
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145 range. The relations between the high and low energy spectral components can probe both 

146 the emission mechanism and the physical conditions in the emission region. The shape of 

147 the high energy spectral energy distribution will be crucial to discriminate between hadronic 
14B cascades and leptonic emission. The spectral breaks at high energy will constrain the Lorentz 

149 factor of the emitting region. Previously undetected emission components might be present 

150 in the light curves such as thermal emission. Finally, temporal analysis of the high energy 

151 delayed component will clarify the nature of the flares seen in the X-ray afterglows. 



2.2.1. Leptonic vs. Hadronic Emission Models 
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It is very probable that particles are accelerated to very high energies close to the emis- 
sion site in GRBs. This could either be in shock fronts, where the Fermi mechanism or 
other plasma instabilities can act, or in magnetic reconnection sites. Two major classes of 
models — synchrotron and inverse Compton emission by relativistic electrons and protons, 
and hadronic cascades — have been proposed for the conversion of particle energy into ob- 
served photon radiation. 

In the leptonic models, synchrotron emission by relativistic ele ctrons can explain the 
10 keV-1 MeV spectrum in ~2/3 of bursts (e.g., see lPreece et al.lll998l ). and inverse Compton 
(IC) scattering of low energy seed photons generally results in GeV band e mission. These pro- 
cesse s could operate in both internal and external shock regions (see, e.g.. lZhang fc Meszaros 
200ll ). with the relativistic electrons in oiie region scattering the 'soft' photons from another 
regio n ( Fragile et al.ll2004 : Fan et aPboOS ; Meszaros et al. 1994 : Waxman 1997 ; Panaitescu et al, 
19981 ). Correlated high and low energy emission is expected if the same electrons radiate 
synchrotron photons and IC scatter soft photons. In Synchrotron Self-Compton (SSC) mod- 
els the electrons' synchrotron photons are the soft photons and thus the high and low en- 
ergy comp onent s should have correlated variability (iGuetta fc Granotll2003l : iGalli &: Guetta 
20081 ). However, SSC models tend to generate a broad uF^, peak in the MeV band, and for 
bursts observed by CG RO this breadth has difficulty accommodating the observed spectra 
( iBaring fc Brabyll2004l ). Fermi, with its broad spectral coverage enabled by the GBM and 
the LAT, is ideally suited for probing this issue further. 



173 
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(Rees & Meszaros 


2005; 


Rvde 


2004, 


2005) 



a component is expected when the outflow becomes optically thin, and would explain low 



energy spectra th a t are too hard for conventional synchrotron models (ICrider et al.l Il997 



Preece et al.lll998l . |2002| ). An additional power law component might underlie this thermal 



component and extend to high energy; this component might be synchrotron emission or 
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179 IC scattering of the thermal photons by relativistic electrons. Fits of the sum of thermal 

180 and power law models to BATSE spectra have been successful (IRy del 12004 |2005| ). but joint 

181 fits of spectra from the two ty pes of GBM detectors a nd the LAT should resolve whether a 

182 thermal component is present ( Battelino et al. 2007a| jbl). 



In hadronic models relativistic protons scatter inelastically off the ~100 keV burst pho- 
tons (j97 interactions) producing (among other possible products) high-energy, neutral pions 
(tt^) that decay, resulting in gamma rays and electrons that then radiate additional gamma 
rays. Similarly, if neutrons in the outflow decouple from protons, inelastic col lisions between 
neutrons and protons can produc e pions and subsequent high energy emission (jPerishev et al. 



2000l : iBahcall fc Meszarod 120001). High energy neutrino s that may be observable are also 



emitted in these interactions (IWaxman fc Bahcalll 119971 ). Many variants of hadronic cas- 
cade models have been proposed: high ene rgy emission from proton- neutron inelastic col- 
lisions early in the evolution of the fireball (IBahcall fc Meszarod 120001): proton-synchrotron 



and photo-meson cascade emission in internal sh ocks (e.g 



Totani 



1998 



Zhang fc Meszaros 



200ll : [Fragile et al.l |2004J: iGupta &: Zhanej 120071 ) ; and proton synchrotron emission in ex- 
ternal shocks^lBo^che^^^ermer 1998 ). A hadronic mode l has been invoked to e xplain 



the additional spectral component observed in ORB 941017 (jPermer fc Atovarul2004f). The 
emiss i on in these models is p redicted to peak in the MeV to GeV band (iBottcher fc Dermer 



19981 : iGupta &: ZhangI 120071 ). and thus would produce a clear signal in the LAT's energy 



band. However, photon -meson interactions would result from a radiatively inefficient fireball 
( IGupta fc ZhangI 120071 ). which is in contrast with the high radiative efficiency that is sug- 
gested by Swift observations (jNousek et al.l l2006l : iGranot et al.ll2006l ). Thus, the hadronic 
mechanisms for gamma-ray production are many, but the Fermi measurements of the tem- 
poral evolution of the highest energy photons will provide strong constraints on these models, 
and moreover discern the existence or otherwise of distinct GeV-band components. 



2.2.2. High-Energy Absorption 



At high energies the outflow itself can become optically thick to photon-photon pair 
production, causing a break in the spectrum. Signatures of internal absorption will constrain 
the bulk Lorentz factor and adiabatic /radiative behavior of the GRB blast wave as a function 



208 of time (iBaring fc Hardindll997l : lLithwick fc Sarill200ll : iGuetta fc Granotll2003l : lBarindl2006 



Granot et al.ll2008l ). Since the outflow might not be steady and may evolve during a burst. 



the breaks should be time- variable, a distinctive property of internal attenuation. Moreover, 
if the attenuated photons and their hard X-ray/soft gamma-ray target photons originate 
from proximate regions in the bursts, the turnovers will approximate broken power-laws. 
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Interestingly, the LAT has aheady provided palpable new advances in terms of constraining 
bulk motion in bursts. For GRB 080916C, the absence of observable attenuation turnovers 
up to around 13 G eV suggests that the bulk Lorentz factor may be well in excess of 500-800 
dAbdo et al-lboooh . 



Spectral cutoffs produced by internal absorption must be distinguished observationally 
from cutoffs caused by interactions with the extragalactic background. The optical depth 
of the Universe to high-energy gamma rays resulting from pair production on infrared and 
optical diffuse extragalactic background radiation can be considerable, thereby preventing 
the radiation from reaching us. These intervening background fields necessarily generate 
quasi-exponential turnovers familiar to TeV blazar studies, which may well be discernible 
from those resulting from internal absorption. Furthermore, their turnover energies should 
not vary with time throughout the burst, another distinction between the two origins for pair 
attenuation. In addition, the turnover energy for external absorption is expected above a few 
lO's of GeV while for internal absorption it may be as low as < 1 GeV ( iGranot et al.ll2008l ). 
Although the external absorption may complicate the study of internal absorption, studies of 
the cutoff as a function of redshift can measur e the universe's optical energy emission out to 
the Population III epoch (with redshift > 7) (Ide Jager fc Steckeii l2002l : ICoppi fc Aharonian 



19971 : lKashhnskvll2005l : iBromm fc Loebll2006f ). 



2.2.3. Delayed GeV Emission 



232 The observations of GRB 940217 (IHurley et al.l Il994j ) demonstrated the existence of 

233 GeV-band emission long after the ~100 keV 'prompt' phase in at least some bursts. With 

234 the multiplicity of shock fronts and with synchrotron and IC components emitted at each 

235 front, many models for this lingering high energy emission are possible. In combination with 

236 the prompt emission observations and afterglow observations by Swift and ground-based 

237 telescopes, the LAT observations may detect spectral and temporal signatures to distinguish 
23B between the different models. 

239 T hese models include: Synchrotron Self-Compton (SSC) emission in late internal shock s 

240 (LIS) Jzhang fc MeszarosI I2OO2I : IWang et al.l bood : iFan et al.l I2OO8I : iGalli fc Gue"ttal l2008h : 

241 external I C (EIC) scattering of LIS photons by the forward shock electrons tha t radiate the 



242 afterglow (IWang et al. 



2001; Granot fc Guetta 



20061) : IC emission in the ex ternal reverse shock (RS) (jWang et al. 



2OO3I: iKobavashi et 



shocks (Meszaros fc Rees 



GaUi fc Piroll2007f ). 



199 



I 



Dermer et al 



20071): and SSC emission in forward external 



2OOOI : IZhang fc MeszarosI I2OOII : IPermeii 12007 
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A high energy IC component may be delayed and have broader time structures relative 
to lower energy components because the scattering ni ay occur in a different region from 
where the soft photons are emitted (IWang et al.l |2006| ) . The correlation of GeV emission 
with X-ray afterg l ow flares observed by Swift would be a diagnostic for different models 
(IWang et allbood : ICaUi fc PirolboO?! : iGaUi fc Guettalboosh . 



2.3. Timing Analysis 
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The LAT's low deadtime and large effective area will permit a detailed study of the 
high energy GRB light curve, which was impossible with the EGRET data as a result of the 
large deadtime that was comparable to typical widths of the peaks in the lightcurve. These 
measures are clearly important for determining the emission region size and the Lorentz 
factor in the emitting fireball. 

The lightcurves of GRBs are frequently complex and diverse. Individual pulses display 
a hard-to-soft evolution, with Ep decreasing exponentially with the burst flux. One method 



of classifying bursts is to examine the spect ral lag, which relates to the de 



high energy and low energy photons (e.g., iNorris et al 



2000 



Folev et al 



Hafizi fc Mochkovitch 



lag v alue indicates hard-to-soft evolution ( Kocevski &: Liangf 12003 
20071 ). i.e., high energy emission arrives earlier than low energy emission. This lag is a direct 



ay in the arrival of 



20081). A positive 



consequence of the spectral evolution of the burst as Ep decays with time. The distributions 
of spectral lags of short and long GRBs a re noticeably different, with the lags of short GRBs 
concentrated in the range ± 30ms (e.g., Norris fc Bonnell 2006 : Yi et al. 2006h. while long 



GRB s have lags coverin g a wid e range with a typical value of 100ms (e.g., iHakkila et al. 
20071 ). IStamatikos et al.l (l2008bl ) study the spectral lags in the Swift data. 



An anti-correlation has bee n discovered betw een the lag and the peak luminosity of 
the GRB at energies ~ 100 keV (INorris et al.ll2000l ). using six BATSE bursts with definitive 
redshift. Brighter long GRBs tend to have a high peak luminosity and short lag, while weaker 
GRBs tend to have lower luminosities and longer lags. This "lag-luminosity relation" has 
been confirmed by using a num ber of Swift GKB s with known redshift (e.g., GRB 060218, 
with a lag greater than 100 s, iLiang et al.l 120061 ). Fermi will be able to determine if this 
relation extends to MeV-GeV energies. 



A subpopulation of local, faint, long lag GRBs has been proposed by iNorrid (120021 ) from 
a study of BATSE bursts, which implies that events with low ' peak fluxes ( FpfSO— 300 keV) ~ 



0.25 ph cm ^ s ^) should be predominantly long lag GRBs. iNorrid (|2002| ) successfully tested 



a prediction that these long lag events are relatively nearby and show some spatial anisotropy. 
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279 and found a concentration towards the lo cal supergalactic plane. This has been confirmed 

280 with the GRBs observed by INTEGRAL JPolev et aliboosj ) where it was found that > 90% 

281 of the weak GRBs with a lag > 0.75 s were concentrated in the supergalactic plan^. 

282 Fermi measures of long lag GRBs will confirm this hypothesis. An underluminous abun- 

283 dant population is infe rred from observations of nearby bursts associated with supernovae 

284 feoderberg et allbooeh . 



285 Moreover, some Quanturn Gravity (QG) theories predict an energy dependent speed-of- 

286 light (see e.g.. iMattinglyl 120051 ) . which is often parameterized as 



v = c{l-{E{z)/E,,)) 



(2) 



287 where E{z) is the photon energy at a given redshift, E{z) = -Eobs(l + z), and Egg is the QG 

288 scale, which may be of order ~ 10^^ GeV. This energy-dependence can be measured from 

289 the difference in the arrival times of different-energy photons that were emitted at the same 

290 time; measurements thus far give E ng greater than a few tirnes 10 ^^ GeV. Such photons might 



291 be emitted in sharp burs t pulses (lAnielino-Camelia et al.lll998l ): measurements have been 

292 attempted (ISchaefeiill999l : iBoggs et al.ll2004l ). The most difficult roadblock to rehable quan- 

293 tum gravity detections or upper limits results from the difficulty in discriminating against 

294 time delays inherent in the emission at the site of the GRB itself, and known to exist from 

295 previous observations. This problem can be addressed by studying a sample of bursts at 

296 different redshifts, or otherwise calibrating this effect. 

297 With the energy difference between the GBM's low energy end and the LAT's high 

298 energy end, the good event timing by both the GBM and the LAT, and the LAT's sensitivity 

299 to high energy photons, the Fermi mission will place interesting limits on Egg. 



possible counterargument has been recently claimed by 



Xiao fc Schaefei 



( 20091) 
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300 3. Description of the Fermi Mission 

301 3.1. Mission Overview 

302 Fermi was launched on June 11, 2008, into a 96.5 min circular orbit 565 km above 

303 the Earth with an inclination of 25.6° to the Earth's equator. During the South Atlantic 

304 Anomaly passages (approximately 17% of the time, on average) the Fermi detectors do not 

305 take scientific data. In FermVs default observing mode the LAT's pointing is offset 35° from 

306 the zenith direction perpendicular to the orbital plane; the pointing will be rocked from one 

307 side of the orbital plane to the other once per orbit. This observing pattern results in fairly 
30B uniform LAT sky exposure over two orbits; the uniformity is increased by the 54 d precession 

309 of the orbital plane. 

310 The mission's telemetry is downlinked 6-8 times per day on the Ku band through the 

311 Tracking and Data Relay Satellite System (TDRSS)o The time between these downlinks, 

312 the transmission time through TDRSS and the processing at the LAT Instrument Science 

313 and Operations Center (LISOC) result in a latency of 6 hours between an observation and 

314 the availability of the resulting LAT data for astrophysical analysis. In addition, when burst 

315 detection software for either detector triggers, messages are sent to the ground through 

316 TDRSS with a ~15 s latency. The mission's burst operations are described in greater detail 

317 below. 



318 3.2. The Large Area Telescope (LAT) 

319 A product of an international collaboration between NASA, DOE and many scientific 

320 institutions across France, Italy, Japan and Sweden, the LAT is a pair conversion telescope 

321 designed to cover the energy band from 20 MeV to greater than 300 GeV. The LAT is 



322 described in greater depth in lAtwood et al.l (120091 ) and here we summarize salient features 



323 useful for understanding the detector's burst capabilities. The LAT consists of an array of 

324 4 x 4 modules, each including a tracker-converter based on Silicon Strip Detector (SSD) 

325 technology and a 8.5 radiation lengths Csl hodoscopic calorimeter. High energy incoming 

326 gamma-rays convert into electron-positron pairs in one of the tungsten layers that are inter- 

327 leaved with the SSD planes; the pairs are then tracked to point back to the original photons' 

328 direction and their energy is measured by the calorimeter. A segmented anti-coincident 

329 shield surrounding the whole detector ensures the necessary background rejection power 



■^See http://msl.jpl.nasa.gov/Programs/tdrss.html 
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330 against charged particles, whose flux outnumbers that of gamma-rays by several orders of 

331 magnitude, and reduce the data volume to fit in the telemetry bandwidth. 

332 Key points of the LAT design are: wide Field-Of-View (FOV — more than 2 sr), large 

333 effective area and excellent Point Spread Function (PSF — see Fig. [2]), short dead time (~ 25 

334 yus per event) and good energy resolution (of the order of 10% in the central region of the 

335 active energy range). As a result, the LAT is the most sensitive high energy gamma-ray 

336 detector ever flown. The study of gamma-ray bursts (GRBs) will take particular advantage 

337 of the improvement in angular resolution — we estimate that two or three photons above 1 
33B GeV will localize a bursts to ~ 5 arcminutes. The reduced dead time will all ow the study 



339 of th e sub-structure of the GRB pulses, typically of the order of milliseconds (IWalker et al. 



20001 ) . with a time resolution that has never before been accessible at GeV energies. 



341 The data telemetered to the ground consists of the signals from different parts of the 

342 LAT; from these signals the ground software must 'reconstruct' the events and filter out 

343 events that are unlikely to be gamma-rays. Therefore, the Instrument Response Functions 

344 (IRFs) depend not only on the hardware but also on the reconstruction and event selection 

345 software. For the same set of reconstructed events trade-offs in the event selection between 

346 retaining gamma rays and rejecting background result in different event classes. There are 

347 currently three standard event classes — the transient, source and diffuse event classes — that 
34B are appropriate for different scientific analyses (as their names suggest). Less severe cuts 

349 increase the photon signal (and hence the effective area) at the expense of an increase in the 

350 non-photon background and a degradation of the PSF and the energy resolution. 

351 The least restrictive class, the transient event class, is designed for bright, transitory 

352 sources that are not background-limited. We expect that the on-ground event rate over the 

353 whole FOV above 100 MeV will be 2 Hz for the transient class and 0.4 Hz for the source 

354 class. In both cases we expect about one non-burst event per minute within the area of the 

355 PSF around the burst position. Consequently, there should be essentially no background 

356 during the prompt emission (with a typical duration of less than a minute) so that the 

357 transient class is the most appropriate — and in fact is the one used for producing all the 
35B results presented in this paper. On the other hand, the analysis of afterglows, which may 

359 linger for a few hours, will need to account for the non-burst background, at least in the low 

360 region of the energy spectrum, where the PSF is larger (see Fig. [2]). 

361 The onboard flight software also performs event reconstructions for the burst trigger. 

362 Because of the available computer resources, the onboard event selection is not as discrim- 

363 inating as the on-ground event selection, and therefore the onboard burst trigger is not as 

364 sensitive because the astrophysical photons are diluted by a larger background flux. Simi- 

365 larly, larger localization uncertainties result from the larger onboard PSF, as shown by the 
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366 left-hand panel of Fig. [2l 

367 3.3. Fermi Gamma-ray Burst Monitor (GBM) 

36B The GBM detects and localizes bursts, and extends FermVs burst spectral sensitivity 

369 to the energy range between 8 keV and 30 MeV or more. It consists of 12 Nal(Tl) (8- 

370 1000 keV) and 2 BGO (0.15-> 30 MeV) crystals read by photomultipliers, arrayed with 

371 different orientations around the spacecraft. The GBM monitors more than 8 sr of the sky, 

372 including the LAT's FOV, and localizes bursts with an accuracy of < 15° {la) onboard, 

373 (< 3° on ground), by comparing the rates in different detectors. The GBM is described in 

374 greater detail in Meegan et al. (2009, submitted). 

375 3.4. FerrnVs Burst Operations 

376 Both the GBM and the LAT have burst triggers. When either instrument triggers, a no- 

377 tice is sent to the ground through the TDRSS within ~ 15 s after the burst was detected and 
37B then disseminated by the Gamma-ray burst Coordinates Network (GCN)0 to observatories 

379 around the world. This initial notice is followed by messages with localizations calculated 

380 by the flight software of each detector. Additional data (e.g., burst and background rates) 

381 are also sent down by the GBM through TDRSS for an improved rapid localization on the 

382 ground by a dedicated processor. 

383 Updated positions are calculated from the full datasets from each detector that are 

384 downlinked with a latency of a few hours. Scientists from both instrument teams analyze 

385 these data, and if warranted by the results, confer. Conclusions from these analyses are 

386 disseminated through GCN Circulars, free-format text that is e-mailed to scientists who 

387 have subscribed to this service. Both Notices and Circulars are posted on the GCN website. 

388 If the observed burst fluxes in either detector exceed pre-set thresholds (which are higher 

389 for bursts detected by the GBM outside the LAT's FOV), the FSW sends a request that the 

390 spacecraft slew to point the LAT at the burst location for a foUowup pointed observation; 

391 currently a 5 hr observation is planned. 

392 In addition to the search for GRB onboard the LAT and manual follow-up analysis by 

393 duty scientists, there is also automated processing of the full science data. This processing 



''Sec http://gcn.gsfc.nasa.gov/ 
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394 performs an independent search for transient events in the LAT data, to greater sensitivity 

395 than is possible onboard, and also performs a counterpart search for all GRB detected within 

396 the LAT FoV. This is described in greater detail in § 16.31 
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Fig. 2. — Left: Comparison of the estimated Point Spread Function (PSF) for the onboard 
and on-ground event reconstruction and selection. The black solid curve is the 68% con- 
tainment angle on-axis for the transient event class, while the dashed curve represents the 
performance of the onboard reconstruction. Right: Comparison of the estimated onboard 
(dashed) and on-ground (solid black curve) on-axis effective areas. These estimates of the 
instrument response are based on simulations of the LAT. 
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Burst Simulations 
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We test the Fermi burst detection and analysis software with simulated data. These 
simulated data are based on our expectations for burst emission in the LAT and GBM 
spectral bands (see § [2]), and on models of the instrument response of these two detectors. 
Since bursts undoubtedly differ from our theoretical expectations, our calculations are more 
reliable in showing the mission's sensitivity to specific bursts than in estimating the number 
of bursts that will be detected. 

We have two 'GR B simulators' that model the burst flux incident on each detector 
( iBattelino et al.l l2007al ). The primary is the phenomenological simulator — described in 
greater detail below in § 14.11 — that dra ws burst parameters from observed distributions. We 



have also created a physical simulator (lOmodei 



2005: Omodei fc Norris 2007: Omodei et al. 



20071 ) that calculates the synchrot ron eraission from the collision of shells in a relativistic 
outflow (the internal shock model — |Piranlll999l ). For a given analysis we assemble an ensem- 
ble of simulated bursts using one of these GRB simulators. To simulate a LAT observation 
of each burst in this ensemble we create a realization of the photon flux, resulting in a 
list of simulated photons incident on the LAT. The LAT's response to this photon flux is 
processed in one of two software paths. The first uses 'GLEAM', which performs a Monte 
Carlo simulation of the propagati on of the photon and its resulting particle shower in the 



LAT (using the GEANT4 too lkit (lAgostinelh et a. 



the different LAT components (lAtwood et al.ll2004l : iBaldini et al.ll2006l ). The photon is then 



2003h) and the det ection of particles in 



'reconstructed' from this simulated instrument response by the same software that processes 
real data. Thus GLEAM maps the incident photons into observed events. Our second, faster, 
processing pathway uses the instrument response functions to map the photons into events 
directly. We note that both approaches use the same input — a list of incident photons-and 
result in the same output — a list of 'observed' events in one of the event classes. In both 
approaches GRBs can be combined with other source types (such as stationary and flar- 
ing AGN, solar flares, supernova remnants, pulsars) to build a very complex model of the 
gamma-ray sky. 

The GRB simulators also provide the input to the GBM simulation software. In this case 
the GRB simulators produce a tim e series of spectral parameters (usually the parameters for 
the 'Band' function — iBandl l2003l — discussed above in § 12. ip . The GBM simulation software 
samples the burst spectrum to create a list of incident photons and then uses a model of 
the GBM response to determine whether each photon is 'detected,' and if so, in which 
energy channel (simulating the GBM's finite spectral resolution). Based on a model from 
the BATSE observations, background counts are added to the burst counts. The GBM 
simulation software outputs count lists, response matrices and background spectra in the 
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433 standard FITS formats used by software such as XSPECO 

434 Because the GRB simulators provide input to both LAT and GEM simulations, sim- 

435 ulated LAT and GBM data can be produced for the same bursts, allowing joint analyses. 

436 The Fermi mission developed the 'Standard Analysis Environment' (SAE) to analyze both 

437 LAT and GBM data. Data can be binned in time, resulting in light curves (see, for example, 
43B Fig. [3]), or in spectra that can be analyzed using a tool such as XSPEC. As will be described 

439 in § [71 joint fits of GBM and LAT data may cover an energy band larger than seven orders of 

440 magnitude (see Fig. [T]). Consequently, Fermi will be a very powerful tool for understanding 

441 the correlation between low-energy and high-energy GRB spectra. 



4.1. Phenomenological Burst Model 



443 
444 
445 
446 



The phenomenological GRB simulator that is used for most of our simulations draws 
from observed spectral and temporal distributions to construct model gamma-ray bursts. 
This modeling assumes that bursts consist of a series of pulses that can be described by a 
universal family of functions (jNorris et al. 1996 ) 



I{t) = A 



exp[-(|t - to\/(Tr) 

exp[-(|t - to\/(Jd) 



t < t 







(3) 



t > to 



447 
44B 
449 



where ar and ad parameterize the rise and decay timescale, and u provides the 'peakiness' of 
the pulse. Although empirically ar ~ 0.33 cr°'*^, we approximate this relation as a^. ~ ad/S. 
The pulse Full Width at Half Maximum (FWHM) is 



W = {ar + ad) ln(2) 



l/u 



(4) 



Pulses are observed to narrow at higher ene rgy in the BATSE energy band (jPavis et al. 



1994J : iNorris et al.lll996l : iFenimore et al.lll995[ ). Although the statistics in the EGRET data 



452 were insufficient to determine whether this narrowing continues in the >100 MeV band, 

453 our phenomenological model assumes that it does. We assume that the FWHM energy 

454 dependence is W{E) oc where ^ is ~0.4 (IFenimore et al.lll995l : iNorris et al.lll996l ). Thus, 

455 we give the pulse shape in eq. [3] an energy dependence by setting 



ad{E) = 0.75 X ln(2)-i/^lVo(E/20 keV)"? 
ar{E) = 0.25 X \n{2)-^/^Wo{E/20 keV)-«. 



(5) 



^See http:/ /heasarc. nasa.gov/xanadu/xspec/ 
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Fig. 3. — Simulated count rate light curve for a BGO detector, two Nal detectors, and the 
LAT for one simulated burst. In this model of the burst spectral evolution, the LAT detects 
counts at the beginning of each pulse; the correlation of the LAT and GBM light curves will 
be a powerful diagnostic of the emission processes. The simulation predicts that the LAT 
would detect a total of 42 gamma rays above 30 MeV in this moderately bright burst of 1 s 
peak flux of 63.37 ph cm~^ s~^ between 30 and 500 keV. 
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456 where Wq is the FWHM at 20 keV. Burst sp ectra in the 10-1000 keV band are well-described 



457 by the 'Band' function (IBand et al.l Il993l ) parameterized in eq. [T] Empirically the Band 

45B function is an adequate description of burst spectra accumulated on short timescales (e.g., 

459 shorter than a pulse width) and over an entire burst. This may be due in part to the poor 

460 spectral resolution of scintillation detectors (such as BATSE and the GBM), but we will treat 

461 this as a physical characteristic of gamma-ray bursts. In the resulting model, the flux /(t, E) 

462 is a product of a Band function with spectral indices a' and (3' and the energy- dependent 

463 pulse shape I(t,E) (eq. [3] with eq. E]) 

f{t,E) = I{t,E) NB.^a{E\No,E„a',(3') ph cm'^ keV~i. (6) 

464 Note that this spectrum is not strictly a Band function because the pulse shape function 

465 does not have a power law energy dependence. 

466 The spectrum integrated over the entire burst is a Band function that is proportional to 

467 the product H^(i?)A^Band(-E|A^O! Ep, a', [3'). Because W{E) is a power law with spectral index 
46B the spectral indices a and (3 for the integrated spectrum are different from the indices for 

469 the instantaneous flux (eq. [6]) 

/oo 
/(t, E)dt = NB.nd{E\No, Ep, a, (3)T = Ao NB.nd{E\No, Ep, a', (3') W{E) 

(7) 

= AoWoNB-,nd{E\No,Ep,a'-^,(3'-0 

470 where T is the burst duration and all the normalizing factors resulting from the integration 

471 are incorporated in Aq. Thus the flux for a single GRB is the sum of many pulses of the 

472 form 

fit, E) = lit, E)iVBand(i? |iVo, Ep,a + i,P + i). (8) 

473 Drawn from observed burst distributions, the same spectral parameters Ep, a and (3 are used 

474 for a given simulated burst. The number of pulses and parameters of each pulse (amplitu de. 



width and peakedness) are also sampled from observed distributions (INorris et al.lll996l ) 



Alternative spectral models have also been simulated; for example, iBattelino et al. 



477 ( l2007al ) describe simulations with a strong thermal photospheric component. 



- 23 - 



47B 5. Semi- Analytical Sensitivity Estimates 

479 The design of the LAT detector provides an ultimate burst sensitivity, regardless of 

480 whether the detection and analysis software achieves this ultimate limit. Thus in this section 

481 we estimate the LAT's burst detection and localization capabilities, and the expected flux 

482 distribution. The following section describes the current burst detection algorithms. 



5.1. Semi-Analytical Estimation of the Burst Detection Sensitivity 



495 
496 
497 



500 
501 



In this subsection we compute the LAT's burst detection sensit ivity using a semi- 
analyt ical approach based on the likelihood ratio test introduced by iNeyman &: Pearson 



486 ( 119281 ). This test is applied extensively to photon-counting experiments ( Cash 



487 has been used to analyze the gamma-ray data from COS-B fjPoUock et al. 



1981 



1979 ) and 



19851 ) and 



488 EGRET (IMattox et al.lll996l ). The statistic for this test is the likelihood for the null hypoth- 

489 esis for the data divided by the likelihood for the alternative hypothesis, here that burst flux 

490 is present. This methodology is the basis of the likelihood tool that will be used to analyze 

491 LAT observations; here we perform a semi-analytic calculation for the simple case of a point 

492 source on a uniform background. 

493 In photon-counting experiments, the natural logarithm of the likelihood for a given 

494 model can be written as 



ln(L) = ln(M() — Aprcd + constant 

photons 



(9) 



where Mj is the predicted photon density at the position and time of ith observed count, 
and Apred is the predicted total number of counts. We compare the log likelihood for the null 
hypothesis that only background counts are present versus the hypothesis that both burst 
and background counts are present. 



The expected number of counts from a burst flux S{E) is 

Ns = Tobs / / A,ff{E)S{E)F{E, dEdQ 

JaQ J El 



(10) 



while the expected number of counts from a background flux B{E) (assumed to be uniformly 
distributed over the sky) is 



pE2 

Nb = Tobs / A,jf{E)B{E)dEAn 
Jei 



(11) 
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502 where A^ff is the effective area and F{E,Q) is the normahzed PSF (which therefore does 

503 not show up in eg. [TTl) . Note that B{E) varies significantly over the sky, but our assumption 

504 is that it is constant over AQ. 

505 The logarithm of the likelihood of the null hypothesis is 

In(Lo) = Tobs [ P A,ff{E) [S{E)F{E, Q) + B{E)] x 
J An J El 

\n{Aeff{E)B{E))dEdVt- Nb • (12) 

506 The actual count rate is assumed to result from both background and burst flux while the 

507 predicted count rates (the Mi in eq. [9] and the total number of counts iVpred) are calculated 
50B only for the background flux (the null hypothesis). 

509 Similarly, the logarithm of the likelihood of the hypothesis that a burst is present is 

'■E2 



ln(Li 



Tabs / / A,ff{E) [S{E)F{E, n) + B{E)] x 
J An J El 



\n{Aeff{E)[S{E)F{E,Q) + B{E)])dEdQ]-{Ns + NB) . (13) 

510 Here both the actual and predicted count rates are calculated for both burst and background 

511 fluxes. 



512 Wilks' theorem (IWilkslll938l ) deflnes the Test Statistic as Ts = — 2(ln(Lo) — In(Li)), and 

513 states that Ts is distributed (asymptotically) as a distribution of m degrees of freedom, 

514 where m is the number of burst parameters. From eqs. [12] and [13] Ts is 

Ts = 2 Tobs [ [ ' Aeff{E)B{E) [(1 + G{E, Q)) In (1 + G{E, Q)) - G{E, Q)] dEdQ (14) 
J An J El 

515 where we have deflned a signal-to-noise ratio G{E,Q) = S{E)F{E,Q)/B{E). 

516 The significance of a source detection in standard deviation units is calculated as N^r = 

517 \^Ts in the case m = 1 {x^ with 1 dof). Here we assume that Wilks' theorem holds, which 
51B might be not absolutely true in a low-count regime (see, in particular, the discussion in 

519 § 16.51) . However, we will see that this method gives a robust estimate of the LAT sensitivity 

520 to GRBs. We can use this method to estimate the LAT sensitivity to GRB. 

521 In our modeling we assume the burst has a 'Band' function spectrum (see eq. [T|) and that 

522 the fiux is constant over a duration Tgrb- Since we seek the optimal detection sensitivity, 

523 we calculate Ts for Tobs = Tgrb- We assume a spatially uniform background with a power 

524 law spectrum 

B{E) = Bo ( ^QQ^^y ) ph cm-2 MeV-^ s'^ si-' (15) 
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525 where the value of the normahzation constant Bq is set to mimic the expected background 

526 rate. For modehng the onboard trigger the background rate above 100 MeV is set to 120 Hz, 

527 while, for the on-ground trigger the background is set to 2 Hz, as will be discussed below. The 
52B spectral index is set to be 7 = —2.1. The results depend on the value of the spectral index; 

529 a detailed study of the dependence of the results as a function of the shape of the residual 

530 background is outside the illustrative goal of this section, thus we omit such discussion. We 

531 require Ts > 25 and at least 10 source counts in the LAT detector, corresponding to a 

532 threshold significance of 5a and a minimum number of GRB counts to see a clear excess 

533 in the LAT data even in the case of very few background events. We use the "transient" 

534 event class described in § 13.21 and compute the minimum 50-300 keV fluence of bursts at 

535 this detection threshold. The burst fluxes in the LAT band depend only on the high energy 

536 power law component of the 'Band' spectrum; assumed values of the low energy power law 

537 spectral index a = —1 and Ep = 500 keV are used to express the spectrum's normalization 
53B in familiar fluence units. Results are shown in Fig. IH at short durations the threshold is 

539 determined by the finite number of burst photons, while the background determines the 

540 threshold for longer durations. This figure predicts that unless other high-energy spectral 

541 components are presen t, the bursts detected by the LAT will be 'hard' with photon indices 



542 /3 near -2 f!Bandll2007l ) 



These estimates consider the detectability of individual bursts. We can compute the 
sensitivity of the LAT detector to GRB considering as input the observ ed distribution of G RB 



545 with known spectral parameters. We use the catalog of bright bursts (iKaneko et al.ll2006l ) to 

546 quantify the characteristics of GRBs. This catalog contains 350 bright GRBs over the entire 

547 life of the BATSE experiment selected for their energy fluence (requiring that the fluence in 
54B the 20-2000 keV band is greater than 2x 10~^ erg/cm^) or on their peak photon flux (over 256 

549 ms, in the 50-300 keV, greater than 10 ph/cm^/s). This subset of burst of the whole BATSE 

550 catalog represents the most comprehensive study of spectral properties of GRB prompt 

551 emission to date and is available electronically from the High-Energy Astrophysics Science 

552 Archive Research Center (HEASARcfl. We restrict our sample of GRB to the ones with a 

553 well reconstructed Epeak] furthermore, we exclude the bursts described by the Comptonized 

554 model (COMP) for which an emission at LAT energy is very unlikely; we also reject bursts 

555 with spectra described by a single power law with undetermined Ep^ak (probably outside the 

556 BATSE energy range). 

557 Considering the field of view of the BATSE experiment and these selection criteria, we 

558 estimate a rate of 50 GRB per year (full sky). For each burst we simulate, the duration, the 



'http://heasarc.gsfc.nasa.gov/ 
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559 energy fluence and the spectral parameters are in agreement with one of the bursts in the 

560 Bright BATSE catalog. Its direction is randomly chosen in the sky, and for each burst we 

561 compute the LAT response functions for that particular direction. Finally, we compute 

562 using eq. [HI The resulting distributions are given by Fig. O 

563 The onboard analysis' larger effective area (Fig. [2]) results in a larger cumulative burst 

564 rate, but not a larger detected rate because of the larger background rate. Events that are 

565 processed onboard by the GRB search algorithm are downloaded, and a looser set of cuts 

566 can be chosen on-ground in order to optimize the signal/noise ratio. We emphasize that this 

567 calculation makes a number of simplifying assumptions. The LAT spectrum is assumed to 
56B be a simple extrapolation of the spectrum observed by BATSE. Spectral evolution within 

569 a burst is not considered. The BATSE burst population was biased by that instrument's 

570 detection characteristics. Nonetheless we estimate that the LAT can detect around 1 burst 

571 per month, with a few bursts per year having more than 100 counts. These few bright bursts 

572 are likely to have a large impact on burst science since detailed spectral analysis will be 

573 possible. 

574 In the framework described in this section, we can also estimate the localization accuracy 

575 for the burst sample, for both onboard and on-ground triggers. If (Jj is the 68% containment 

576 radius for the single photon PSF, then the localization is computed as 



^GRB 



577 that, in terms of the previously defined quantities, is 



_ ITgrb A,„(E)S(E) 



578 The factor of 3 takes into account the non-gaussianity of the PSF, and was estimated by 



579 iBurnettI (120071 ). We compute the localization accuracy for each burst in our sample. Fig. [6] 

580 shows the results. In each plot the detected burst are represented by red triangles, while the 

581 blue empty circles are the bursts with LAT counts that did not pass our detection condition. 

582 These results show that the LAT can localize bursts with sub-degree accuracy, both 

583 onboard and on-ground. The GRB yield is greater and bursts are better-localized on-ground 

584 than onboard. The on-ground analysis is available only after the full dataset is downlinked 

585 and processed. This process can lasts few hours, depending on the position of the downlink 

586 contact. Onboard localization is delivered quasi- real time with onboard alerts. For those 

587 bursts, multiwavelength follow-ups will be feasible for bursts localized within a few tens of 

588 arcminutes. For example, the FOV of Swift's XRT is about 0.4° and is of the same order 
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Fig. 4. — Threshold fluence as a function of the GRB duration, for on-ground detection and 
for on-axis incidence. Threshold fluence increases by factor of ~ 2 for z-axis angles of 50 
degrees. Different lines are related to different spectral index. Also plotted are the observed 
bursts from the BATSE catalog. 




Fig. 5. — Integrated number of GRBs per year as a function of the number of LAT counts. 
The solid curve shows all bursts in the sample, while the dashed curve gives the detected 
bursts. Left panel: on-ground analysis ("transient" class, 2 Hz background rate above 
100 MeV). Right panel: onboard analysis (120 Hz background rate). 
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589 as the FOV of the typical mid-size optical or near-IR (NIR) telescope. Afterglow searches 

590 in the optical and NIR are very successful — ~60% of the Swift bursts have been associated 

591 with optical and NIR afterglows. Fig. E] shows that a sizeable fraction of Fermi GRB 

592 detections will be localized within these requirements, and relatively large FOV ground- 

593 based observatories (~30 arcmin) with optical/NIR filters (I, z, J, H, K) should produce a 

594 fairly high detection rate for the afterglows of LAT-detected GRBs. 
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Fig. 6. — Number of LAT counts vs. localization accuracy. In each panel the red triangles 
denote detected bursts and the open blue circles show undetected bursts. The left and 
right panels are for the on-ground and onboard localizations. Thus the on-ground analysis 
results in a slightly larger burst detection rate and a better localizations. The superior track 
reconstruction and background reduction outweighs the smaller effective area in increasing 
the on-ground detection rate. 
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5.2. Estimated LAT Flux Distribution 



596 We now consider the full GRB model described in § IHfor estimating the expected LAT 

597 flux distribution. This is, of course, very dependent on the assumptions of the GRB model, 

598 and the final result should be considered only as a prediction of the flux distribution. 
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We use the bright BATSE catalog ( IKaneko et al.l l2006l ) for the burst population, as 
described in the previous section. In addition, we also select a sub-sample of bursts for 
which beta is more negative than -2. This is motivated by the fact that a power law index 
greater than -2 implies a divergence in the released content of energy, thus those value are 
unphysical and a cut-off should take place. The measurements yielding beta greater than -2 
are questionable and suggest either an ill-determined quantity for a true spectrum that is 
in reality softer, or an additional spectral break above the energies measured with BATSE. 
Given the duration, the number of pulses is flxed by the total burst duration. Pulses are 
combined together in order to obtain a flnal Tgo duration. Correlations between duration, 
intensity, and spectral parameters are automatically taken into account as each of these 
bursts corresponds to an entry in the Kaneko et al. catalog. The emission is extended up to 
high energy with the model described in § |H 

We emphasize again that this model ignores possible intrinsic cutoffs (resulting from 
the high end of the particle distribution or internal opacity — § 12.2.21) . and additional high- 
energy components suggested by the EGRET observations (§ 12.11) . High-energy emission 
(>10 GeV) is also sensitive to cosmological attenuation due to pair production between the 
GRB radiation and the Extragalactic Background Light (EBL — § 12.2.21) . The uncertain EBL 
spectral energy distribution resulting from the absence of high redshift data provides a variety 
of theoretical models for such diffuse radiation. Thus the observation of the high-energy cut- 
off as a function of the GRB distance can, in principle, constr ain the background l ight. In 
our simulation we include this effect, adopting the EBL model in lKneiske et al.l (120041 ). Short 
bursts are thought to be the result of the merging of compact objects i n bin ary systems, 
so we adopt the short burst redshift distribution from iGuetta fc PiranI (120051 ) , while long 
bursts are related to the explosi ve end of massive stars, whose distributions are well traced 
by the Star Formation History (jPorciani &: Madaull200ll ). 



In Fig. [7] the sampled distributions are shown. The Dashed line histogram is obtained 
from the full bright burst BATSE catalog. In order to increase the number of burst in 
the fleld of view of the LAT detector we over-sampled the original catalog by a factor 1.4. 
The dark flUed histograms show the distribution of GRB with at least 1 count in the LAT 
detector, and the light filled histograms are the sub-sample of detected GRB with beta < -2. 



We simulate approximately ten years of observations in scanning mode. The orbit of 
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Fig. 7. — Parameter distributions for the simulated bursts of the bright burst BATSE catalog 
(dashed lines). Filled dark histograms represent the GRBs with more than 1 predicted count 
above 100 MeV in the LAT detector, while for the light filled histograms we have also required 
that the high-energy spectral index beta is more negative than -2. The distributions show 
the logarithm of the duration, the fluence, the peak flux distribution, the low and high energy 
spectral indexes and the logarithm of the energy of the peak of the uFi, spectrum. 
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the Fermi satellite, the South Atlantic Anomaly (SAA) passages and Earth occultations 
are all considered. In Fig. [8] we plot the number of expected bursts per year as a function 
of the number of photons per burst detected by the LAT. The different couples of lines 
refer to different energy thresholds (100 MeV, 1 GeV, and 10 GeV). Dashed lines are the 
same computation but using only the sub-sample of GRBs with beta more negative than 
-2 (the light filled distribution in Fig. [7]). The EBL attenuation affects only the high- 
energy curve, as expected from the theory, leaving the sensitivities almost unchanged below 
10 GeV. Assuming that the emission component observed in the 10-1000 MeV band continues 
unbroken into the LAT energy band, we estimate that the LAT will independently detect 
approximately 10 bursts per year, depending on the sensitivity of the detection algorithm; 
approximately one burst every three months will have more than a hundred counts in the 
LAT detector above 100 MeV: these are the bursts for which a detailed spectral or even time 
resolved spectral analysis will be possible. If we restrict our analysis to the sub-sample of 
bursts with beta more negative than -2, these numbers decrease. Nevertheless, even if we 
adopt this conservative approach, LAT should be able to detect independently approximately 
1 burst every two months, and will be able to detect radiation up to tens of GeV. 

With the assumed high-energy emission model a few bursts per year will show high- 
energy prompt emission, with photons above 10 GeV. These rates are in agree ment with the 



numb er of bursts dete cted in the LAT da ta after few mo nths (GRB0808 25C (iBouvier et al. 



20081), GRB080916C flTaiima et all 120081 ). GRB081024B flOmodeill2008[ )). but the statistics 
is still low for any strong constraint on the burst population. 
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Number Of Photons Detected 

Fig. 8. — Model-dependent LAT GRB sensitivity. The GRB spectrum is extrapolated from 
BATSE to LAT energies. The all-sky burst rate is assumed to be 50 GRB yr~-^ full sky 
(above the peak flux in 256 ms of 10 ph s~^ cm^^ in the 50-300 keV or with an energy flux 
in the 20-2000 keV band greater than 2x 10^^ erg/cm^), based on BATSE catalog of bright 
bursts. The effect of the EBL absorption is included. Different curves refer to different 
energy thresholds. Dashed curves are the result of the analysis excluding very hard bursts, 
with a beta greater than -2. 
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6. Gamma-Ray Burst Detection 



652 The rapid detection and localization of bursts is a major goal of the Fermi mission. 

653 Both Fermi instruments will search for bursts both onboard and on-ground. These searches 

654 will detect bursts on different timescales and with different sensitivities. Here we focus on 

655 LAT burst detection, but for completeness we describe briefly GBM burst detection. 

656 6.1. GBM Burst Detection 

657 Onboard the Fermi observatory the GBM will use rate triggers that monitor the count 

658 rate from each detector for a statistically significant increase. Similar to the BATSE detec- 

659 tors, the GBM as a whole will trigger when two or more detectors trigger. A rate trigger 

660 compares the number of counts in an energy band AE over a time bin At to the expected 

661 number of background counts in this AE-At bin; the background is estimated from the rate 

662 before the time bin being tested. The GBM trigger uses the twelve Nal detectors with vari- 

663 ous energy bands, including A£'=50-300 keV, and time bins from 16 ms to 16.384 s. Note 

664 that the BATSE trigger had one energy band — usually Ai?=50-300 keV — and the three time 

665 bins At =0.064, 0.256, and 1.024 s. The GBM burst detection algorithms are described in 

666 greater detail in Meegan et al. (2009, submitted). 

667 When the GBM triggers it sends a series of burst alert packets through the spacecraft and 

668 TDRSS to the Earth. Some of these burst packets, including the burst location calculated 

669 onboard, will also be sent to the LAT to assist in the LAT's onboard burst detection. Burst 

670 locations are calculated by comparing the rates in the different detectors; each the detectors' 

671 effective area varies across the FOV. In addition, the GBM will send a signal over a dedicated 

672 cable to the LAT; this signal will only inform the LAT that the GBM has triggered. 

673 The continuous GBM data that are routinely telemetered to the ground can also be 

674 searched for bursts that did not trigger the GBM onboard. These data will provide rates 

675 for all the GBM detectors in 8 energy channels with 0.256 s resolution and in 128 energy 

676 channels with 4.096 s resolution. In particular, if a burst triggers the LAT but not the GBM, 

677 these rates will at the very least provide upper limits on the burst flux in the GBM energy 

678 band. 
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679 6.2. Onboard LAT Detection 

680 The LAT flight software will detect bursts, localize them, and report their positions 

681 to the ground through the burst alert telemetry. The rapid notification of ground-based 

682 telescopes through GCN will result in multi-wavelength afterglow observati ons of GRBs 



683 with k nown high energy emission. The onboard burst trigger is described in iKuehn et al. 



(120071). 



685 The onboard processing that results in the detection of a GRB can be subdivided into 

686 three steps: initial event filtering; event track reconstruction; and finally burst detection and 

687 localization. In the first step all events — photons and charged particles — that trigger the 

688 LAT hardware are filtered to remove events that are of no further scientific interest. The 

689 events that survive this first filtering constitute the science data stream that is downlinked 

690 to the ground for further processing. These events are also fed into the second step of the 

691 onboard burst processing pathway. 

692 The second step of the burst pathway attempts to reconstruct tracks for all the events in 

693 the science data stream using the 'hits' in the tracker's silicon strip detectors that indicate the 

694 passage of a charged particle. The burst trigger algorithm uses both spatial and temporal 

695 information, and therefore a 3-dimensional track that points back to a photon's origin is 

696 required. Tracks can be calculated for only about a third the events that are input to this 

697 step, although surprisingly the onboard track-finding efficiency is 80% to 90% of the more 

698 sophisticated ground calculation. However, the onboard reconstruction is less accurate, 

699 resulting in a larger PSF onboard than on-ground, as is shown by Fig. [2l A larger fraction of 

700 the incident photons survive the onboard filtering than survive the on-ground processing at 

701 the expense of a much higher non-photon background onboard than on-ground; consequently 

702 the onboard effective area is actually larger than the on-ground effective Fig. [2] shows. 

703 The rate of events that pass the onboard gamma filter (currently the same event set 

704 that is downlinked and thus available on-ground) is ~400 Hz. The rate that events are 

705 sent to the onboard burst trigger, which requires 3-dimensional tracks, is ~120 Hz. The 

706 on-ground processing creates a transient event class with a rate of ~2 Hz. Thus onboard the 

707 burst trigger must find a burst signal against a background of ~120 non-burst events, while 

708 on-ground this background is only ~2 Hz. This difference in non-burst background rate sets 

709 fundamental limits on the onboard and on ground burst detection sensitivities. 

710 The third step in the burst processing is burst detection, which considers the events 

711 that have passed all the filters of the first two steps, and thus have arrival times, energies 

712 and origins on the sky. When a detector such as the GBM provides only event rates, the 

713 burst trigger can only be based on a statistically significant increase in these rates. However, 
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714 when a detector such as the LAT provides both spatial and temporal information for each 

715 event, then an efficient burst trigger will search for temporal and spatial event clustering. 

716 Most searches for transients bin the events in time and space (if relevant) , but the LAT uses 

717 an unbinned method. 

718 The LAT burst trigger searches for statistically significant clusters in time and space. 

719 The trigger has two tiers. The first tier identifies potentially interesting event clusters for 

720 further investigation by the second tier; the threshold for the first tier allows many false 

721 tier 1 triggers that arc then rejected by the second tier. The first tier operates continuously, 

722 except while the second tier code is running. A GBM trigger is equivalent to a first tier 

723 trigger in that the GBM's trigger time and position arc passed directly to the second tier. 

724 Tier 1 operates on sets of N events that survived the first two steps, where currently 

725 is in the range of 40-200. The effective time window that is searched is A^ divided by the 

726 event rate; for an event rate of 120 Hz and these values of A^, the time window is 1/3-5/3 s. 

727 Each of these A^ events is considered as the seed for a cluster consisting of all events that 

728 are within 9q of the seed; currently Oq — 17°, approximately the 68% containment radius of 

729 the onboard 3D tracks at low event energies. A clustering statistic, described below, is then 

730 calculated for each cluster. A tier 1 trigger results when a clustering statistic for any cluster 

731 exceeds a threshold value. A candidate burst location is then calculated from the events of 

732 the cluster that resulted in the tier 1 trigger. 

733 The onboard burst localization algorithm uses a weighted average of the positions of the 

734 cluster's events. The weighting is the inverse of the angular distance of an event from the 

735 burst position. Since the purpose of the algorithm is to find the burst position, the averaging 

736 must be iterated, with the weighting used in one step calculated from the position from the 

737 previous step. The initial location is the unweighted average of the events positions. The 

738 convergence criterion is a change of 1 arcmin between iterations (with a maximum of 10 

739 iterations). The position uncertainty depends on the number and energies of events, but the 

740 goal is an uncertainty less that 1°. Using Monte Carlo simulations, this methodology was 

741 found to be superior to others that were tried. 

742 The tier 1 trigger time and localization (or if the GBM triggered, its trigger time and 

743 burst position) are then passed to the second tier. Because the second tier is run relatively 

744 infrequently, it can consider a much larger set of events than the first tier. Currently 500 

745 events are considered, which corresponds to a time window of ~4.2 s. A cluster is then 

746 formed from all events in this set that are within 02 (~ 10°) of the tier 1 burst location. A 

747 clustering statistic is then calculated for this cluster, and if its value exceeds a threshold, a 

748 tier 2 trigger results and the cluster events are run through the localization algorithm. The 

749 resulting trigger time, burst location and number of counts in four energy bands are then 
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750 sent to the ground through the burst alert telemetry. The second tier is run repeatedly after 

751 a tier 1 trigger in case the burst brightens resulting in a larger cluster centered on the tier 1 

752 position, and consequently a tier 2 trigger (if one has not yet occurred) and a better burst 

753 localization (if a tier 2 trigger does occur) . 

754 The clustering statistic is based on the probabilities that the cluster's events have the 

755 observed distances from the chistcr seed position and the arrival time separations, under the 

756 null hypothesis that a burst is not occurring. Assuming events are thrown uniformly onto a 

757 sphere (the null hypothesis), the probability Ps of finding an event within 6 degrees of the 

758 cluster seed position is 

1 — cos(^) , , 

Ps = 18 

1 - cos(6l^) 

759 where it is assumed that there are no events at more than 9m — 115° (the performance is 

760 not sensitive to this parameter). Thus for a cluster of M events the spatial contribution to 

761 the clustering statistic is 



(19) 



762 The temporal part of the cluster probability assumes that the event arrival time follows 

763 a Poisson distribution (again the null hypothesis). The probability that the arrival times of 

764 two subsequent events differ by AT is 

pt = l-exp[-rtAT] , (20) 

765 where Tt is the rate at which events occur within the area of the cluster. The temporal 

766 contribution of each cluster to the clustering statistic is 

M M 



Pt^YI llogiobtJI = E Fogio(l - e"'*"""^') ■ (21) 



i=l 



767 The trigger criterion is 

^PT + Ps>e (22) 

768 where ^ is an adjustable parameter that assigns relative weights to the spatial and temporal 

769 clustering, and © is the threshold. The two tiers may use different values of both ^ and ©. 

770 The overall false trigger rate depends on the tier 2 value of ©. 

771 The parameters used by the onboard burst detection and localization software are sen- 

772 sitive to the actual event rates, and will ultimately be set based on flight experience. Cur- 

773 rently the thresholds are set high enough to preclude any triggers, and diagnostic data is 
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774 being downlinked and studied. The thresholds will eventually be lowered, keeping the false 

775 trigger rate at an acceptable level. 

776 Based on preliminary calculations using a burst population based on BATSE, we es- 

777 timate ~1 bursts every two months will be detected and localized to 1° (see Fig. [5] and 
77B Fig.E]). 
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779 6.3. LAT Ground-Based Blind Search 

780 A burst detection algoritlim will be applied on the ground to all LAT counts after 

781 the events are reconstructed and classified to detect bursts that were not detected by the 

782 onboard algorithm, the GBM, or other missions and telescopes. Thus this 'blind search' is 

783 similar to the first tier of the onboard burst detection algorithm. The ground-based search 

784 will be performed after each satellite downlink; to capture bursts that straddle the downlink 

785 boundaries, some counts from the previous downlink are buffered and used in searching for 

786 bursts in the data from a given downlink. The ground-based blind search algorithm is very 

787 similar to the onboard algorithm described in the previous section, but will benefit from the 

788 full ground-based event reconstruction and background rejection techniques that are applied 

789 to produce the LAT counts used for astrophysical analysis. For these data, the particle 

790 background rates will be lower than the onboard rates by at least two orders-of-magnitude. 

791 Furthermore, the reconstructed photon directions and energies will be more accurate than 

792 the onboard quantities. Fig. [2] compares the 68% containment angle as a function of the 

793 photon energy for the onboard and on-ground LAT count datasets. 

794 In addition to differing in the reconstruction and background filtering, the ground- 

795 based analysis treats the input data slightly differently. The first stage of the ground- 

796 based algorithm is apphed to consecutive sets of 20 to 100 counts. As with the onboard 

797 algorithm, the number of counts analyzed is configurable and will be adjusted with the 

798 growth of our knowledge of GRB prompt emission in the LAT band and of the residual 

799 instrumental background. However, in contrast to the onboard algorithm, the data sets do 

800 not overlap. This ensures that each segment is statistically independent and generally better 

801 separates the log-probability distributions of the null case (i.e., where there is no burst) 

802 from the distributions computed when burst photons are present. Fig. [H] shows the reference 

803 distribution for the null case derived from simulated background data. We modeled the low 

804 end (large negative values) of the distribution with a Gaussian, and set the burst detection 

805 threshold at 5a from the fitted peak. Since this distribution is derived from pre-launch Monte 

806 Carlo simulations with assumed incident particle distributions and other expected on-orbit 

807 conditions, the thresholds are being re-calibrated with real fiight data. Since we perform 

808 an empirical threshold calibration, we can neglect the constant normalization factors in the 

809 denominators of the single event probabilities shown in eqs. [18] and [201 

810 The overall log-probability is the sum of spatial and temporal components (see eq. [22]) . 

811 which we weight equally (^=1). Fig. ]T0] shows the 2D distributions for the temporal and 

812 spatial components. The dashed line in Fig. [10] corresponds to the 5a threshold with this 

813 weighting. Fig. [TT] shows the time history of the log-probabilities as applied to the GRB grid 

814 data. The excursions across the threshold line indicate the burst candidates. 
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log-probability 

Fig. 9. — Distribution of log-probability values under the null hypothesis obtained from 
applying the ground-based version of the GRB search algorithm to sets of 20 counts. The 
shaded region indicates the range over which a Gaussian function, shown in red, was fit to 
these data. The resulting 5a threshold at an overall log-probability value of —117 is plotted 
as the vertical dashed line. Burst candidates are required to have log-probabilities below 
this threshold. 
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Fig. 10. — 2D histogram of the spatial and temporal log-probability components. The dashed 
line indicates the 5—0" threshold (an overall log-probability value of —117) derived from the 
null distribution (figure [9]). Burst candidates are required to lie below this line. 
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Fig. 11. — Time history of the ground-based log-probabihty. The horizontal dashed hne 
shows the Sa threshold derived from the Gaussian function fit to the log-probabilities dis- 
tribution under the null hypothesis (Fig. [9]). Burst candidates are required to lie below this 
line. 
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B15 While the onboard burst trigger performs two passes through the data with the temporal- 

816 spatial clustering likelihood algorithm, the ground-based detection analysis performs only one 

817 such pass. If a candidate burst is found in the ground-based analysis, counts from a time 

818 range bracketing the trigger time undergoes further processing to determine the significance 

819 of the burst. If the burst is sufficiently significant, it is localized and its spectrum is analyzed. 

820 These analyses use the unbinned maximum likelihood method that is applied to LAT point 

821 sources. 



822 6.4. GRB Candidate Follow-up Processing 

823 When a candidate burst location and trigger time is provided by the ground-based 

824 blind search, a LAT or GBM onboard trigger, or another burst detector such as Swift — 

825 we will call this a first stage detection — a LAT ISOC data processing pipeline will analyze 

826 the LAT counts to determine the significance of a possible LAT detection. This step in 

827 deciding whether the LAT has detected a burst is similar to the tier two analysis of the 

828 onboard algorithm. If the LAT has detected a burst, the pipeline will localize the burst and 

829 determine its temporal start and stop. All of the analyses described in this section will be 

830 performed using the "transient" class. These data selections have a larger effective area at 

831 a cost of somewhat higher instrumental background, particularly in the 50-200 MeV range. 

832 For bright transients, such as are expected for GRBs, this trade-off is advantageous given 

833 the short time scales. 

834 The first step in the follow-up processing is determining the time interval straddling the 

835 candidate burst during which the LAT count rate is greater than the expected background 

836 rate. The counts are selected from a 15° acceptance cone centered on the candidate burst 

837 position and from a 200 second time window centered on the candidate burst trigger time. 

838 This time window is designed to capture possible precursor emission that may be present in 

839 the LAT band. Both the acceptance cone radius and the time window size are configurable 

840 parameters in the processing pipeline. With this acceptance cone radius, the total event rate 

841 from non-GRB sources is expected to be < 0.1 to 0.5 Hz for normal scanning observations, 

842 depending on how far the candidate position is from the brightest parts of the Galactic 

843 plane emission. Th e event arriva l times are analyzed using a Bayesian Blocks algorithm 



844 ( iJackson et al.l l2003l : IScargld Il998l ) that aggregates arrival times in blocks of constant rate 

845 and identifies "change points" between blocks with statistically significant changes in event 

846 rate. The burst start and stop time are identified as the first and last change points from 

847 the resulting light curve. An example of the results of this analysis is shown in Fig. [T2l 



848 



If no change points are found within the 200 second bracketing time window, then the 
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B49 counts from the first stage time window and burst position will be used in calculating upper 

850 limits. In these cases, the position refinement step will be skipped and background model 

851 components will be included in the significance and upper limits analysis. 

852 If application of the Bayesian Block algorithm to the LAT arrival times finds a statisti- 

853 cally significant increase in the count rate above background, i.e., if at least two change points 

854 were found, then further analysis uses only the counts between the first and last change 

855 points to exclude background. The position is refined with the standard LAT maximum 

856 likelihood software that folds a parameterized input source model through the instrument 

857 response functions to obtain a predicted distribution of observed counts. The parameters of 
B5B the source model are adjusted to maximize the log-likelihood of the data given the model. 

859 For these data, the background counts are sufficiently small that a model with the different 

860 background components usually used in point source analysis is not needed, and a model with 

861 a single point source should suffice to localize the burst. The burst spectral parameters and 

862 burst coordinates are adjusted within the extraction region to maximize the log-likelihood, 

863 and the best-fit position is thereby obtained. Error contours are derived by mapping the 

864 likelihood surface in position space, with 90% confidence limit (CL) uncertainties given by 

865 the contour corresponding to a change in the log-likelihood of 2.305. This value is equal to 

866 for 2 degrees-of-freedom (dof). Fig. [13] shows an example counts map with the 90% 

867 CL contour overlaid. 

868 For spectral analysis and the definitive burst significance calculation we use the counts 

869 within the first and last change points and at the center of a 15° radius acceptance cone 

870 around the maximum likelihood position. Again we use maximum likelihood to derive the 

871 basic burst parameters from the LAT data alone. Since this is an automated procedure, a 

872 simple power-law model is chosen as the default. For brighter bursts, background model 
B73 components are not needed. For fainter bursts, such as those burst candidates for which 

874 we only have a first stage detection, including the background is essential to determine the 

875 significance of a faint burst in the LAT data and for deriving upper limits. 

876 6.5. Quantifying Significance and Upper Limits 

877 As discussed in § 15. H the likelihood ratio test (LRT) is a natural framework for hypoth- 

878 esis testing, and we will use this method for quantifying the significance of a candidate burst. 

879 The background models used for the null hypothesis (i.e., that a burst is not present) can 

880 be simplified considering the expected number of counts from each background component 

881 over the short GRB time scales (< (9(10^) s). For determining the significance of a source, 

882 we compute the test statistic defined in eq. [TH We are fairly conservative and require a 
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Fig. 12. — LAT counts light curve for a simulated burst (solid histogram) and a piece- wise 
constant light curve derived using the Bayesian Blocks analysis of the event arrival times 
(dashed histogram). 




Fig. 13. — LAT counts map for the simulated burst in Fig. [12] using only the counts between 
the first and last change points. The best-fit position and 90% error contour derived from 
the maximum likelihood analysis are overlaid. The color scale on the right shows the counts 
per pixel. 
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B83 Ts > 25, corresponding to 5cr for 1 dof, in order to claim a detection. 

884 Upper limits may be computed in several different ways. A method that has been used 

885 in the past for GRBs and other transient astronomical sources is a variant of the classical 

886 "on source-off source" measurement. In this method, one defines an appropriate background 

887 interval prior to the time of the candidate burst, and using the inferred background levels, 

888 one derives an upper-limit for the source flux given the counts that are observed during 

889 the interval containing the candidate burst. Application of this procedure requires that 

890 the observing conditions (instrument response, intrinsic background rates, etc.) during the 

891 background interval be sufficiently similar to those for the interval containing the putative 

892 signal. For the short time spans appropriate for GRBs (< 100 s), simulations have shown 

893 that the instrumental background rates are fairly constant; in survey mode, at fixed rocking 

894 angle, the LAT FOV scans across the sky at a few degrees per minute, so the instrument 

895 response to a given source location will be roughly constant as well. A major benefit of 

896 this procedure is that it is model- independent. However, being model independent, it is also 

897 fairly conservative; and in general, it will not give the most constraining upper-limit. 

898 A more stringent upper-limit may be computed with the "profile likelihood" method. 

899 In this method the normalization of the source flux (or a parameter that determines this 

900 normalization) is varied while fitting all the other model parameters, resulting in the variation 

901 of the log-likelihood (the fitting statistic) as a function of the source normalization. For a 

902 two-sided interval, under Wilks' theorem the 90% confidence region corresponds to a change 

903 in the log-likelihood from the extremum of 2.71/2, i.e., = Ax^/2 for 1 dof. For a one-sided 

904 interval, as in the case of an upper-limit, this corresponds to a 95% CL. 

905 To illustrate the method, we apply this analysis to simulated data. Fig. [T4l shows a LAT 

906 counts map and lightcurve for the time and location of a simulated burst that was detected 

907 in the GEM, but is not evident in the LAT data. The best-fit flux and error estimate for a 

908 point source is 3.2 ±4.5 x 10~^ph cm~^ s^^ for energies E > 100 MeV. The test statistic for 

909 the point source is Tg = 0.67, consistent with the flux measurement's large error bars and 

910 the lack of a burst detection. Fig. [15] shows the fitted counts spectrum and residuals from 

911 this fit. Fig. [TBI shows the change in log-likelihood as a function of scanned fiux value. For 

912 a 95% CL upper limit, we find a value of 1.3 x 10"^ ph cm~^ s~^. 

913 To check the method's validity, we ran Monte Carlo simulations under the same ob- 

914 serving conditions and using the source model and best-fit parameters from the likelihood 

915 analysis as inputs, and we analyzed each simulation to find the best-fit fiux. The left panel 

916 of Fig. [T7] shows the distribution of fitted fiuxes for these simulations, and the right panel 

917 shows the normalized cumulative distribution for these data and the cumulative distribution 

918 inferred by computing the corresponding probability from the profile likelihood curve 
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919 shown in Fig. [161 




Fig. 14. — Left: LAT counts map for a 60 s time window containing the GBM trigger time 
of a simulated burst. The GBM location and 4.5° error circle are plotted. The dashed line 
indicates the location of the Galactic plane. The color scale on the right shows the counts 
per pixel. Right: Counts light curve for these data. The GBM trigger time is indicated by 
the vertical dashed line. 
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Fig. 15. — Fitted counts spectrum and residuals for the data shown in Fig. [TH The contri- 
butions of the three model components are plotted as the long dashed curves, and from top 
to bottom, are the Galactic diffuse, extragalactic diffuse, and point source. The solid curve 
is the sum of the three components. 



-51 - 




0.4 0.6 0.8 1.0 1.2 1.4 

Flux (ph/cm^2/s) 

x10 ' 



Fig. 16. — Change in the log-hkehhood as a function of GRB flux for E > 100 MeV. The 
horizontal dashed hne indicate the 95% CL corresponding to an upper-hmit of 1.3 x 10~^ 
ph cm~^ s~^. 
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Fig. 17. — Left: Distribution of fitted fluxes for the point source representing the GRB 
derived from 766 LAT simulations using the best-flt model obtained from the original dataset. 
Right: The solid curve is the normalized cumulative distribution determined from the fltted 
flux distribution. The dotted curve is the cumulative fraction that would be predicted by 
the likelihood proflle shown in flgure [161 
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7. Spectral Analysis 



921 To demonstrate the spectral analysis that will be possible with the Fermi data, we 

922 present two sample analyses, the first the joint fit of GBM and LAT count spectra, and the 

923 second the search for a cutoff in the LAT energy band. In both cases we use transient class 

924 LAT counts. In general, bursts are short but bright, and thus we can tolerate the higher 

925 background rate of the transient class to increase the number of burst counts. While we focus 

926 here on LAT-GBM joint fits, such fits will also be possible between the Fer mi detectors and 

927 those of other missions, such as Swift (IStamatikos et al.ll2008al : lBandll2008l ). 



7.1. GBM and LAT Combined Analysis 



929 In this example, we assume that a simulated burst was detected and localized by the 

930 GBM. Analysis of the LAT data found 160 transient event class photons in a 20° region 

931 surrounding the GBM position during the 3 s prompt phase observed by the GBM; the 

932 Automated Science Processing (ASP) that will be run after the LAT events are reconstructed 

933 ( §3.2p localized the burst with an uncertainty radius of 0.05°. Fig. [TS] shows the GBM and 

934 LAT light curves. 

935 The simulated GBM and LAT data, both event lists, were accumulated over the burst's 

936 prompt phase, and the LAT events were binned into 10 energy bins. Two Nal and one BGO 

937 detector provided count spectra. The GBM background spectra used to simulate the counts 
93B were used as the background for the GBM count spectra, while the LAT data were assumed 

939 not to be contaminated by background events. We performed a joint fit to the 4 count 

940 spectra (from 2 Nal, one BGO and t he LAT det ectors) with the standard X-ray analysis 

941 tool XPSEC using the Cash statistic (ICashlll979l ). The 'Band' spectrum (eq. [1]) was used 
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Fig. 18. — GBM Nal (left) and LAT (right) light curves of the prompt emission from the 
simulated burst. 
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942 to create the simulated data and for the joint fit. Fig. [19] shows the simulated data (with 

943 error bars) and best-fit model (histogram). The fit yielded a = —0.97 ±0.05 (input value of 

944 -1.09) and (3 = -1.80 ± 0.01 (input value of -1.90). 

945 Thus Fermi will measure the energy spectrum of bursts over 7 orders of magnitude 

946 in energy through its combination of detectors. The energy bands of the Nal and BGO 

947 detectors overlap in the energy region of the peak energy, and the BGO and the LAT energy 
94B bands also overlap. 



949 7.2. Study of GRB high-energy properties with the LAT 

950 Whether the burst spectrum is a simple power law in the LAT energy band, or has a 

951 cutoff spectrum is of great theoretical interest (see § 12.2.21) . Therefore, we simulated and 

952 then fit spectra with such cutoffs to determine if they would be detectable. 

953 We used the simulation software described in § 14.11 to simulate 5 years of Fermi ob- 

954 servations. In this simulation, the temporal and spectral properties of GRBs were based 

955 on a phenomenological or physical model, including not only synchrotron emission but also 

956 inverse Compton emission for a few bursts. The simulated spectra did not have any intrinsic 

957 cutoffs, but included gamma-ray absorption by the Ex tragalactic Backgrou nd Light (EBL) 



95B between the burst and the Earth, following the model of iKneiske et al.l (120041 ). This extrinsic 

959 cut-off only appears at the highest energies (at least 10 GeV), depending on the distance of 

960 the bursts. 

961 The search for high-energy cut-offs was performed using only simulated LAT data. First 

962 we selected those bursts that have no inverse Compton component, and more than 20 LAT 

963 counts. Each count spectrum was fit both by a simple power law and by a power law with 

964 an exponential cutoff with characteristic energy E^. 

965 The likelihoods of the resulting fits were examined to evaluate the improvement of the fit 

966 by adding the cutoff (one additional parameter). The difference of the likelihoods follows a 

967 x^-distribution with one degree of freedom, with the null hypothesis probability distribution 
96B shown in Fig. [20l Two bursts exhibit a very small probability of being consistent with no 

969 cutoff, and thus we consider these bursts to have a statistically significant high-energy cutoff. 

970 While both bursts have average redshifts (1.71 and 3.35) compared to the full sample, they 

971 are very bright, with more than 1000 photons detected. 



972 For these two bursts we per formed a second fit using the parameterisation of the EBL 

973 cut-off proposed by iReyed (120071 ) where the cutoff is exp(— r), with r = 1 + {E — Ei)/P for 
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Fig. 19. — Photon spectrum of the simulated burst: in the top panel, crosses show the data 
of the different sub-detectors (two Nal detectors in black and red, one BGO in green, and 
the LAT in blue) and the histogram denotes the best fit of a Band function. The bottom 
panel shows the ratio of the simulated data to the fit model. 
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Fig. 20. — x^-probability of the difference of the hkehhoods of fits of a power law with and 
without an exponential cutoff: a probability of < 5.7 x 10"'' corresponds to a 5a detection 
of a cutoff. 
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974 E > El — P, and otherwise; Ei is the redshift-dependent energy where the optical depth 

975 is unity, and P is a redshift-dependent energy scahng factor. The two fitted values of Ei 

976 (Sl.Slgg GeV and 43.51'^q [J GcV), arc in good agreement with the true values (46.6 GeV and 

977 30.7 GcV) of the model used for the simulation. Thus the LAT will be sensitive to cutoffs 

978 in the brightest bursts, with good spectral reconstruction. 
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8. Coordination with Other Burst Missions 



8.1. Coordination with Swift 
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The Fermi detectors will provide few localizations accurate to less than 10 arcmin that 
are necessary for t he optical foUowups that can determine redshifts. On the other hand, the 
Swift instruments (iGehrels et al.ll2004l ) — the Burst Alert Telescope (BAT), the X-Ray Tele- 
scope (XRT) and the Ultraviolet-Optical Telescope (UVOT) — provide progressively better 
burst localizations that are rapidly disseminated by the GRB Coordinate Network (GCN), 
resulting in multiwavelength followup observations and frequently burst redshifts. How- 
ever, the BAT's 15-150 keV energy band is often insufficient to determine the spectrum 
of the prompt burst emission, particularly Ep, the 'peak energy' where most of the burst 
energy is radiated (see § 14.11) : Ep is important not only for burst energetics but also for re 



ported relationships between intrinsic burst parameters ( lAmatill2006l : iGhirlanda et al.l 12004 



Firmani et aPbooeh that may turn bursts into standard candles. And for those bursts where 



the BAT can determine Ep, Swift cannot determine whether there is a second emission com- 
ponent above the 15-150 keV band (as discussed in §2.21) . In addition. Swift's burst afterglow 
observations 'only' extend to the X-ray band {E < 10 keV); as discussed in ^[2.11 E GRET 



detected GeV-band prompt and afterglow emissions (IHurlev et al 



1994 



Dingusll2003h . Thus 



Fermi and Swift capabilities complement each other (jstamatikos et aL 2008a ): between the 



UVOT, XRT, BAT, GBM and LAT, the two mission's observations span 11 energy decades. 

The Fermi and Swift missions are working to increase the number of bursts that are ob- 
served simultaneously by the BAT and the LAT; this will increase the number of bursts with 
localizations, redshifts, spectra and optical through gamma-ray afterglows. Simultaneous 
burst observations by Fermi and other burst missions (e.g., AGILE, INTEGRAL, Konus- 
Wind, RHESSI, Suzaku-WAM) will also c omplement each other and permit cross-calibration, 
but Swift's pointing is the most flexible (lBandll2008l ). 



Fermi's and Swift's low earth orbits (altitudes of ~565 and ~590 km, respectively) are 
inclined to the Earth's equator by 25.6° and 20.6°, respectively. The two orbits will beat 
with a period of ~13 days, that is, the two missions will be on the same side, or opposite 
sides, of the Earth with a nearly two week period. Because of the uniformity of the LAT's 
sky-exposure and the large FOVs of the BAT and the LAT, the relative inclination of the 
two orbits (which can be as small as 5° or as large as 46°) has little effect on the overlap of 
the FOVs. The relative inclination varies with a period of approximately 6.5 years. 

In general Fermi will survey the sky, pointing the LAT 35° above or below the orbital 
plane (as described in § |3]). On the other hand, every orbit Swift points the Narrow-Field 
Instruments (NFIs — the XRT and UVOT) at a number of targets that satisfy the mission's 
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1014 observational constraints: the NFIs cannot be pointed near the Sun, moon, horizon or ram di- 

1015 rection; anti-Sun observations are preferred to increase the detection of bursts during Earth's 

1016 night. Since FermVs observing mode will not change, but Swiffs timeline is by design ex- 

1017 tremely flexible, increasing the overlap between the mission's FOVs, and thus increasing the 

1018 number of simultaneous burst detections, will be done through Swift's targeting. Between 

1019 following-up bursts the Swift NFIs are used for other observation programs (and will observe 

1020 Fermi sources). By choosing NFI targets at times that will increase the LAT-BAT overlap, 

1021 we estimate that this overlap can be improved by a factor of ~2 without sacrificing Swiffs 

1022 science objectives. Note that increasing the BAT-LAT overlap will by necessity increase the 

1023 overlap between the BAT and GBM. 

1024 Swift detects ~100 bursts per year, and approximately one LAT detection per month is 

1025 anticipated, although this prediction of the LAT's detection rate is based on extrapolations 

1026 from lower energy (see § 15. 2p . Given the differences in the detectability of typical bursts, we 

1027 assume that Swiffs BAT will detect all the bursts that the LAT will detect when the burst 

1028 is in both their FOVs The LAT's larger FOV compensates for the BAT's greater ability 

1029 to detect typical bursts, resulting in comparable detection rates. Based on a number of 

1030 modeling assumptions, and assuming that Swiffs targeting can increase the overlap of the 

1031 BAT and LAT FOVs by x2, we estimate ~10 BAT bursts per year with LAT detections 

1032 or upper limits, and ~4 LAT bursts per year with BAT detections. We emphasize that our 

1033 estimates of the LAT detection rate assumes that the 10-1000 keV component observed by 

1034 BATSE, BAT and now the GBM extrapolates unbroken into the LAT's energy band. 

1035 8.2. TeV Observations 

1036 The synergy between Fermi and ground-based telescopes operating above a few tens of 

1037 GeV will expand the study of the still-unknown spectral and temporal properties of GRBs 

1038 above a few GeV. Extending the analysis of burst temporal and spectral properties to even 

1039 higher energies would have a large impact on the knowledge of the particle acceleration and 

1040 emission processes occurring in the burst environment. High energy spectra would probe 

1041 the distant Universe, revealing the universe's transparency to high-energy gamma-rays and 

1042 measuring EBL. The requirements for a good coordination of Fermi with TeV observatories 

1043 are quite simple, and we examine the potential of such simultaneous observations in terms 

1044 of expected rates of alerts and sensitivity. 

1045 Major TeV observatories operate above ~100 GeV (or somewhat lower for the next 

1046 generation of instruments), and Imaging Atmospheric Cherenkov Telescopes (lACTs) have 

1047 a sensitivity of 10~^^ to 10~^ erg cm~^ to the latter part of the prompt phase and early 



- 60 - 



1048 afterglow emission of GRBs (i.e., from ~10 s to a few hours after the trigger time). The 

1049 observatories' duty cycle, FOV and sky coverage will determine their response to Fermi 

1050 alerts. With a high duty cycle (~100%) and a good sky coverage (~20%), ground arrays 

1051 like MILAGRO and ARGO will be able to react to any alert provided by the GBM or the 

1052 LAT. In contrast, lACTs like CANGAROO, HESS, MAGIC, VERITAS, or STACEE have 

1053 a low duty cycle (~10%) because they observe only during clear and moonless nights, but 

1054 they can slew to any location within a few minutes and access ~20% of the sky. Because of 

1055 their small FOV (~5°), lACTs will require a GRB position accuracy of ±1° and thus will 

1056 respond effectively to LAT alerts only. 

1057 Using a phenomenological model to describe GRB properties in the LAT range, we 

1058 combine the estimated GRB detection rate (1 GRB per month) with the above duty cycle 

1059 and sky coverage to compute the possible joint observations by Fermi and TeV experiments. 

1060 Fermi should provide ~40 alerts (including 2 to 5 LAT alerts) per year during the prompt 

1061 burst phase, that ground arrays will be able to follow up. Few of them will be followed- up 

1062 by lACTs due to localization accuracy and to observing time constraints. The LAT detected 

1063 bursts per year suitable for TeV foUowup should be considered as the highest priority targets 

1064 in TeV telescope plans. A few afterglows per year may be also followed-up by lACTs, while 

1065 ground arrays will probably be much less sensitive to afterglows. 



8.3. Neutrino Observations 



A major step forward in understanding of the microphysics of the GRB central engines 
might be achieved via the detection of non-electromagnetic emission such as gravitational 
waves (lAbbott et al.ll2005l ) and neutrinos. Because they are weakly-interacting, neutrinos 
are unique (albeit elusive) cosmic messengers because they are not absorbed nor deflected 
on th eir way to the observer. The viability of high energy neutrino astronomy (iGaisser et al. 



19951 ) opens a new observing channel that complements the high energy electromagnetic 



1073 spectrum that will be probed directly by the LAT. 



Hadronic fireball models ( §2.2.ip . predict a taxonomy of correlated MeV to EeV neutri- 
nos of varying flavor and arrival times. Ideal for detection are ~TeV-PeV muon neutrinos 
(IWaxman fc Bahcalll 119971 ) produced as the leptonic decay products of photomeson interac- 

[n] — s> /i+ + z/„ — s> e"*" + i/g + ^'h + ^'u) within the internal shocks 



1077 tions {p + J ^ A" 



TT ' 



of the relativistic fireball. Since the prompt gamma rays act as the ambient photon target 
field, the burst neutrinos are expected to be spatially and temporally coincident with the 
gamma-ray emission. Therefore Antar ctic Cherenkoy teles copes such as Antarctic Muon and 
Neutrino Detector Array (AMANDA) flAhrens et al.ll2002l ) and IceCube (lAhrens et al.ll2004l ) 
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can perform a nearly background-free search for burst neutrinos correl ated with the prompt 
gamma-ray emission (IStamatikos et al.ll2005l : IStamatikos fc Bandll2006l ) . Neutrino telescopes 
have FOVs determined by their position on the Earth, and accumulate and preserve their 
data, and therefore need not to respond to bursts in realtime. Instead, the neutrino data 
archived is searched periodically for neutrinos correlated with the time and position of prompt 
burst emission. Analysis of AMANDA data has resulted in t he most stringent uppe r limits 
upon correlated multi- flavored neutrino emission from GRBs (lAchterberg et al.ll2007l . l2008l ). 
AMANDA'S km-scale successor, IceCube, is currently under construction with anticipated 
completion by ~2010, and thus will operate during the Fermi era. 
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9. Conclusions and Future Work 

In this paper we provided an overview of the LAT's capabilities to reveal the rich burst 
phenomenology in the >100 MeV band at which the EGRET observations merely hinted, 
and which theoretical scenarios predict. These capabilities can be realized only through 
efficient analysis techniques and software. In this final section we discuss the future analysis 
development that we anticipate during the early part of the Fermi mission. 

Burst triggers are applied to the LAT data both onboard and on-ground. The onboard 
trigger contends with a higher non-burst background rate, but can provide burst notifications 
and localizations within tens of seconds after the burst, while the on-ground trigger is more 
sensitive because the background can be reduced, but the burst notification and localizations 
have a ~3 hr latency. The thresholds for both triggers depend on the actual instrument 
response and background rates that are only now being evaluated. Thus during the mission's 
early phase we will tune the detection algorithms to minimize false triggers and maximize 
the detection sensitivity. 

In particular, we are investigating various 'cuts' of the reconstructed events used by the 
on-ground detection algorithms. These cuts do not merely increase or decreased the effective 
area and the background rate, but also change their energy dependence. Relative changes in 
the effective area and background rate affect the detectability of bursts of different durations, 
since the background is less important for detecting short bursts. 

The GBM and LAT spectra will be analyzed jointly, giving spectral fits from ~8 keV to 
over 300 GeV, a bandpass of up to 7.5 energy decades. Typically the spectral analysis will 
fit the parameters of functional forms such as the 'Band' function. 

However, given the theoretical uncertainties in the underlying GRB spectrum in the 
LAT band (e.g., the unknown high energy attenuation by the EBL and intrinsic photon 
fields), we will explore model-independent spectral reconstruction. Deconvolution of instru- 
ment response effects in the Poisson statistic s regime is notorious l y diff icult, but there have 



been advances in recent years. For example, iNowak fc KolaczykI (120001 ) derived a Bayesian 



multiscale framework that is inspired by wavelet methods, but adapted for Poisson statis- 
tics; using these methods , they reconstruct ed a Solar flare emission line spectrum observed 



1124 



by CGRO^s COMPTEL. iD'Agostinil (119951 ) derived another Bayesian iterative method for 
deconvolving spectra; uncertainties on the unfolded distribution can be estimated from a 
covariance matrix. 

Thus we anticipate an exciting mission exploring new burst phenomena and developing 
the techniques to extract the maximum information from the LAT. 
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1128 following the launch of Fermi. His presence on the Fermi team is already greatly missed. 

1129 We thank the members of the LAT instrument team, GBM instrument team and the 

1130 Fermi Project for their exceptional efforts in developing the Fermi observatory. M. Sta- 

1131 matikos is supported by an NPP Fellowship at NASA-GSFC administered by ORAU. 

1132 The Fermi LAT Collaboration acknowledges support from a number of agencies and 

1133 institutes for both the development and the operation of the LAT as well as scientific data 

1134 analysis. These include the National Aeronautics and Space Administration and the De- 
lias partment of Energy in the United States, the Commissariat a I'Energie Atomique and the 

1136 Centre National de la Recherche Scientifique / Institut National de Physique Nucleaire et de 

1137 Physique des Particules in France, the Agenzia Spaziale Italiana and the Istituto Nazionale 

1138 di Fisica Nucleare in Italy, the Ministry of Education, Culture, Sports, Science and Technol- 

1139 ogy (MEXT), High Energy Accelerator Research Organization (KEK) and Japan Aerospace 

1140 Exploration Agency (JAXA) in Japan, and the K. A. Wallenberg Foundation, the Swedish 

1141 Research Council and the Swedish National Space Board in Sweden. Additional support 

1142 from the Istituto Nazionale di Astrofisica in Italy for science analysis during the operations 

1143 phase is also gratefully acknowledged. 

REFERENCES 

1145 Abbott, B., et al. 2005, Phys. Rev. D, 72, 042002 

1146 Abdo, A., et al. 2009, Science, 323, 1688 

1147 Achterberg, A., et al. 2007, Ap. J., 664, 397 

1148 — . 2008, Ap. J., 674, 357 

1149 Agostinelli, S., et al. 2003, Nuclear Instruments and Methods in Physics Research A, 506, 

1150 250 

1151 Ahrens, J., et al. 2002, Phys. Rev. D, 66, 012005 

1152 — . 2004, Astroparticle Physics, 20, 507 

1153 Amati, L. 2006, MNRAS, 372, 233 



-64- 



1154 Amelino-Camelia, G., Ellis, J., Mavromatos, N. E., Nanopoulos, D. V., & Sarkar, S. 1998, 

1155 Nature, 395, 525 

1156 Atwood, W., et al. 2004, in Calorimetry in Particle Physics, ed. C. Cecchi, P. Cenci, P. Lu- 

1157 brano, & M. Pepe, 329-336 

1158 Atwood, W., et al. 2009, Ap. J., 697, 1071 

1159 Bahcall, J. N., & Meszaros, P. 2000, Physical Review Letters, 85, 1362 

1160 Baldini, L., et al. 2006, Nucl. Phys. Proc. Suppl., 150, 62 

1161 Band, D. 2003, Ap.J., 588, 945 

1162 Band, D., Matteson, J., Ford, L., Schaefer, B., Palmer, D., Teegardcn, B., Cline, T., Briggs, 

1163 M., Paciesas, W., Pendleton, G., Fishman, G., Kouveliotou, C., Meegan, C., Wilson, 

1164 R., & Lestrade, P. 1993, Ap. J., 413, 281 

1165 Band, D. L. 2007, in American Institute of Physics Conference Series, Vol. 921, The First 

1166 GLAST Symposium, ed. S. Ritz, P. Michelson, & C. A. Meegan, 446-447 

1167 Band, D. L. 2008, in American Institute of Physics Conference Series, Vol. 1000, American 

1168 Institute of Physics Conference Series, 121-124 

1169 Baring, M. G. 2006, Ap. J., 650, 1004 

1170 Baring, M. G., & Braby, M. L. 2004, Ap. J., 613, 460 

1171 Baring, M. G., & Harding, A. K. 1997, Ap. J. Lett., 481, L85+ 

1172 Battelino, M., Ryde, F., Omodei, N., & Band, D. L. 2007a, in American Institute of Physics 

1173 Conference Series, Vol. 921, American Institute of Physics Conference Series, ed. 

1174 S. Ritz, P. Michelson, & C. A. Meegan, 478-479 

1175 Battelino, M., Ryde, F., Omodei, N., & Longo, F. 2007b, in American Institute of Physics 

1176 Conference Series, Vol. 906, American Institute of Physics Conference Series, ed. 

1177 M. Axelsson & F. Ryde, 28-39 

1178 Boggs, S. E., Wunderer, C. B., Hurley, K., & Coburn, W. 2004, Ap. J. Lett., 611, L77 

1179 Bottcher, M., & Dermer, C. D. 1998, Ap. J. Lett., 499, L131+ 



-65- 



1180 Bouvier, A., Band, D., Bregeon, J., Chiang, J., Cutini, S., Dingus, B., Gehrels, N., Pukazawa, 

1181 Y., Hayashida, M., Longo, F., McEnery, J., Ohno, M., Omodei, N., Pelassa, V., 

1182 Piron, P., Sanchez, D., Scargle, J., Tajima, H., Tanaka, T., & Thayer, G. 2008, GRB 

1183 Coordinates Network, 8183, 1 

1184 Briggs, M. S. 1999, in Astronomical Society of the Pacific Conference Series, Vol. 190, 

1185 Gamma-Ray Bursts: The Pirst Three Minutes, ed. J. Poutanen & R. Svensson, 133 — h 

1186 Bromm, V., & Loeb, A. 2006, Ap. J., 642, 382 

1187 Burnett, T. H. 2007, in American Institute of Physics Conference Series, Vol. 921, American 

1188 Institute of Physics Conference Series, ed. S. Ritz, P. Michelson, & C. A. Meegan, 

1189 530-531 

1190 Cash, W. 1979, Ap. J., 228, 939 

1191 Coppi, P. S., & Aharonian, F. A. 1997, Ap. J. Lett., 487, L9+ 

1192 Crider, A., Liang, E. P., Smith, I. A., Preece, R. D., Briggs, M. S., Pendleton, G. N., 

1193 Paciesas, W. S., Band, D. L., & Matteson, J. L. 1997, Ap. J. Lett., 479, L39+ 

1194 D'Agostini, G. 1995, Nuclear Instruments and Methods in Physics Resarch A, 362, 487 

1195 Davis, S. P., Norris, J. P., Kouveliotou, C, Fishman, G. J., Meegan, C. A., & Paciesas, 

1196 W. S. 1994, in American Institute of Physics Conference Series, Vol. 307, Gamma- 

1197 Ray Bursts, ed. G. J. Fishman, 182 — h 

1198 de Jager, O. C, & Stecker, F. W. 2002, Ap. J., 566, 738 

1199 Derishev, E. V., Kocharovsky, V. V., & Kocharovsky v., V. 2000, in American Institute of 

1200 Physics Conference Series, Vol. 526, Gamma-ray Bursts, 5th Huntsville Symposium, 

1201 ed. R. M. Kippen, R. S. Mallozzi, & G. J. Fishman, 460-464 

1202 Dermer, C. D. 2007, ArXiv Astrophysics e-prints 

1203 Dermer, C. D., & Atoyan, A. 2004, Astron. & Astrophys., 418, L5 

1204 Dermer, C. D., Chiang, J., & Mitman, K. E. 2000, Ap. J., 537, 785 

1205 Dingus, B. L. 2003, in American Institute of Physics Conference Series, Vol. 662, Gamma- 

1206 Ray Burst and Afterglow Astronomy 2001: A Workshop Celebrating the First Year 

1207 of the HETE Mission, ed. G. R. Ricker & R. K. Vanderspek, 240-243 

1208 Fan, Y.-Z., Piran, T., Narayan, R., & Wei, D.-M. 2008, Mon. Not. RAS, 384, 1483 



-66- 



1209 Fan, Y. Z., Zhang, B., & Wei, D. M. 2005, Ap. J., 629, 334 

1210 Fenimore, E. E., in 't Zand, J. J. M., Norris, J. P., Bonnell, J. T., & Nemiroff, R. J. 1995, 

1211 Ap. J. Lett., 448, L101+ 

1212 Firmani, C, Ghisellini, G. andAvila-Reese, V., & Ghirlanda, G. 2006, MNRAS, 370, 185 

1213 Foley, S., McGlynn, S., Hanlon, L., McBreen, S., & McBreen, B. 2008, Astron. & Astrophys., 

1214 484, 143 

1215 Ford, L. A., Band, D. L., Mattcson, J. L., Briggs, M. S., Pendleton, G. N., Preece, R. D., 

1216 Paciesas, W. S., Teegarden, B. J., Palmer, D. M., Schaefer, B. E., Cline, T. L., 

1217 Fishman, G. J., Kouveliotou, C, Meegan, C. A., Wilson, R. B., & Lestrade, J. P. 

1218 1995, Ap. J., 439, 307 

1219 Fragile, P. C, Mathews, G. J., Poirier, J., & Totani, T. 2004, Astroparticle Physics, 20, 591 

1220 Gaisser, T. K., Halzen, F., & Stanev, T. 1995, Physics Reports, 258, 173 

1221 Galh, A., & Guetta, D. 2008, Astron. & Astrophys., 480, 5 

1222 Galh, A., & Piro, L. 2007, Astron. & Astrophys., 475, 421 

1223 Gehrels, N., Chincarini, G., Giommi, P., Mason, K. O., Nousek, J. A., Wells, A. A., White, 



1224 N. E., Barthelmy, S. D., Burrows, D. N., Cominsky, L. R., Hurley, K. C., MarshaU, 

1225 F. E., Meszaros, P., Roming, P. W. A., Angelini, L., Barbier, L. M., Belloni, T., 

1226 Campana, S., Caraveo, P. A., Chester, M. M., Citterio, O., Cline, T. L., Cropper, 

1227 M. S., Cummings, J. R., Dean, A. J., Feigelson, E. D., Fenimore, E. E., Frail, D. A., 

1228 Fruchter, A. S., Garmire, G. P., Gendreau, K., Ghisellini, G., Greiner, J., Hill, J. E., 

1229 Hunsberger, S. D., Krimm, H. A., Kulkarni, S. R., Kumar, P., Lebrun, F., Lloyd- 

1230 Ronning, N. M., Markwardt, C. B., Mattson, B. J., Mushotzky, R. F., Norris, J. P., 

1231 Osborne, J., Paczynski, B., Palmer, D. M., Park, H.-S., Parsons, A. M., Paul, J., 

1232 Rees, M. J., Reynolds, C. S., Rhoads, J. E., Sasseen, T. P., Schaefer, B. E., Short, 



1233 A. T., Smale, A. P., Smith, I. A., Stella, L., Taghaferri, G., Takahashi, T., Tashiro, 

1234 M., Townsley, L. K., Tueller, J., Turner, M. J. L., Vietri, M., Voges, W., Ward, M. J., 

1235 Willingale, R., Zerbi, F. M., & Zhang, W. W. 2004, Ap.J., 611, 1005 

1236 Ghirlanda, G., Ghisellini, G., & Lazzati, D. 2004, Ap. J., 616, 331 

1237 Giuliani, A., Mereghetti, S., Fornari, F., Del Monte, E., Feroci, M., Marisaldi, M., Esposito, 

1238 P., Perotti, F., Tavani, M., Argan, A., Barbiellini, G., Boffelli, F., Bulgarelli, A., 

1239 Caraveo, P., Cattaneo, P. W., Chen, A. W., Costa, E., D'Ammando, F., di Cocco, G., 



-67- 



1240 Donnarumma, I., Evangelista, Y., Fiorini, M., Fuschino, P., Galli, M., Gianotti, F., 

1241 Labanti, C, Lapshov, I., Lazzarotto, P., Lipari, P., Longo, P., Morselli, A., Pacciani, 

1242 L., Pellizzoni, A., Piano, G., Picozza, P., Prest, M., Pucella, G., Rapisarda, M., 

1243 Rappoldi, A., Sofiitta, P., Trifoglio, M., Trois, A., Vallazza, E., Vercellone, S., Zanello, 

1244 D., Salotti, L., Cutini, S., Pittori, C, Preger, B., Santolamazza, P., Verrccchia, P., 

1245 Gehrels, N., Page, K., Burrows, D., Rossi, A., Hurley, K., Mitrofanov, I., & Boynton, 



1246 W. 2008, Astron. & Astrophys., 491, L25 

1247 Gonzalez, M. M., Dingus, B. L., Kaneko, Y., Preece, R. D., Dermer, C. D., & Briggs, M. S. 

1248 2003, Nature, 424, 749 

1249 Granot, J., Cohen-Tanugi, J., & do Couto e Silva, E. 2008, Ap. J., 677, 92 

1250 Granot, J., & Guetta, D. 2003, Ap. J. Lett., 598, 111 

1251 Granot, J., Konigl, A., & Piran, T. 2006, Mon. Not. RAS, 370, 1946 

1252 Guetta, D., & Granot, J. 2003, Ap. J., 585, 885 

1253 Guetta, D., & Piran, T. 2005, Astron. & Astrophys., 435, 421 

1254 Gupta, N., & Zhang, B. 2007, Mon. Not. RAS, 380, 78 

1255 Hafizi, M., & Mochkovitch, R. 2007, Astron. & Astrophys., 465, 67 

1256 Hakkila, J., Gibhn, T. W., Young, K. C, Fuller, S. P., Peters, C. D., Nolan, C, Sonnett, 

1257 S. M., Haglin, D. J., & Roiger, R. J. 2007, Ap. J. Supp., 169, 62 



1258 Hoover, A. S., Kippen, R. M., & McConnell, M. L. 2005, Nuovo Cimento C Geophysics 

1259 Space Physics C, 28, 825 

1260 Hurley, K., Dingus, B. L., Mukhcrjce, R., Sreekumar, P., Kouveliotou, C., Meegan, C., 

1261 Fishman, G. J., Band, D., Ford, L., Bertsch, D., Cline, T., Fichtcl, C., Hartman, R., 

1262 Hunter, S., Thompson, D. J., Kanbach, G., Maycr-Hasselwandcr, H., von Montigny, 

1263 C., Sommer, M., Lin, Y., Nolan, P., Michelson, P., Kniffen, D., Mattox, J., Schneid, 

1264 E., Boer, M., k Niel, M. 1994, Nature, 372, 652 

1265 Jackson, B., Scarglc, J. D., Barnes, D., Arabhi, S., Alt, A., Gioumousis, P., Gwin, E., 

1266 Sangtrakulcharoen, P., Tan, L., & Tsai, T. T. 2003, ArXiv Mathematics e-prints 

1267 Kaneko, Y., Gonzalez, M. M., Preece, R. D., Dingus, B. L., & Briggs, M. S. 2008, Ap. J., 

1268 677, 1168 



-68- 



1269 Kaneko, Y., Preece, R. D., Briggs, M. S., Paciesas, W. S., Meegan, C. A., & Band, D. L. 

1270 2006, Ap. J. Supp., 166, 298 

1271 Kashlinsky, A. 2005, Ap. J. Lett., 633, L5 

1272 Kneiske, T. M., Bretz, T., Mannheim, K., & Hartmann, D. H. 2004, Astron. & Astrophys., 

1273 413, 807 

1274 Kobayashi, S., Zhang, B., Meszaros, P., & Burrows, D. 2007, Ap. J., 655, 391 

1275 Kocevski, D., & Liang, E. 2003, Ap. J., 594, 385 

1276 Kouvehotou, C, Preece, R., Bhat, N., Fishman, G. J., Meegan, C. A., Horack, J. M., 



1277 Briggs, M. S., Paciesas, W. S., Pendleton, G. N., Band, D., Matteson, J., Palmer, D., 

1278 Teegarden, B., & Norris, J. P. 1994, Ap. J. Lett., 422, L59 

1279 Kuehn, F., Bonnell, J., Hughes, R., Norris, J., Ritz, S., Russell, J., Smith, P., & Winer, B. 

1280 2007, in American Institute of Physics Gonference Series, Vol. 921, American Institute 

1281 of Physics Gonference Series, ed. S. Ritz, P. Michelson, & G. A. Meegan, 556-557 



1282 Liang, E.-W., Zhang, B.-B., Stamatikos, M., Zhang, B., Norris, J., Gehrels, N., Zhang, J., 

1283 & Dai, Z. G. 2006, Ap. J. Lett., 653, L81 

1284 Lithwick, Y., & Sari, R. 2001, Ap. J., 555, 540 

1285 Mattingly, D. 2005, Living Reviews in Relativity, 8, 5 

1286 Mattox, J. R., Bertsch, D. L., Ghiang, J., Dingus, B. L., Digel, S. W., Esposito, J. A., 



1287 Fierro, J. M., Hartman, R. G., Hunter, S. D., Kanbach, G., Kniffen, D. A., Lin, 

1288 Y. G., Macomb, D. J., Mayer-Hasselwander, H. A., Michelson, P. F., von Montigny, 

1289 G., Mukherjee, R., Nolan, P. L., Ramanamurthy, P. V., Schneid, E., Sreekumar, P., 

1290 Thompson, D. J., & WiUis, T. D. 1996, Astrophysical Journal v.461, 461, 396 



1291 Meszaros, P., & Rees, M. J. 1994, Mon. Not. RAS, 269, L41+ 

1292 Meszaros, P., Rees, M. J., & Papathanassiou, H. 1994, Ap. J., 432, 181 

1293 Neyman, J., & Pearson, E. 1928, Biometrika 

1294 Norris, J. P. 2002, Ap. J., 579, 386 

1295 Norris, J. P., & Bonnell, J. T. 2006, Ap. J., 643, 266 

1296 Norris, J. P., Marani, G. F., & Bonnell, J. T. 2000, Ap. J., 534, 248 



-69- 



1297 Norris, J. P., Nemiroff, R. J., Bonnell, J. T., Scargle, J. D., Kouveliotou, C, Paciesas, W. S., 

1298 Meegan, C. A., & Fishman, G. J. 1996, Ap. J., 459, 393 

1299 Nousek, J. A., Kouveliotou, C, Grupe, D., Page, K. L., Granot, J., Ramirez- Ruiz, E., Patel, 

1300 S. K., Burrows, D. N., Mangano, V., Barthelmy, S., Beardmore, A. P., Campana, S., 

1301 Capalbi, M., Chincarini, G., Cusumano, G., Falcone, A. D., Gehrels, N., Giommi, P., 

1302 Goad, M. R., Godet, O., Hurkett, C. P., Kennea, J. A., Moretti, A., O'Brien, P. T., 

1303 Osborne, J. P., Romano, P., Tagliaferri, G., & Wells, A. A. 2006, Ap. J., 642, 389 

1304 Nowak, R. D., & Kolaczyk, E. D. 2000, IEEE Transactions on Information Theory, 46, 1811 

1305 Omodei, N. 2005, in High Energy Gamma-ray Experiments, ed. A. De Angelis & O. Mansutti, 

1306 189-196 

1307 Omodei, N. 2008, GRB Coordinates Network, 8407, 1 

1308 Omodei, N., & Norris, J. 2007, in American Institute of Physics Conference Series, Vol. 921, 

1309 American Institute of Physics Conference Series, ed. S. Ritz, P. Michelson, & C. A. 

1310 Meegan, 472-475 

1311 Omodei, N., et al. 2007, AIP Conf. Proc, 906, 1 

1312 Paciesas, W. S., Meegan, C. A., Pendleton, G. N., Briggs, M. S., Kouvehotou, C, Koshut, 

1313 T. M., Lestrade, J. P., McCoUough, M. L., Brainerd, J. J., Hakkila, J., Henze, W., 

1314 Preece, R. D., Connaughton, V., Kippen, R. M., Mallozzi, R. S., Fishman, G. J., 

1315 Richardson, G. A., & Sahi, M. 1999, Ap. J. Supp., 122, 465 

1316 Panaitescu, A., Meszaros, P., & Rees, M. J. 1998, Ap. J., 503, 314 

1317 Piran, T. 1999, Physics Reports, 314, 575 

1318 Pollock, A. M. T., Bennett, K., Bignami, G. F., Bloemen, J. B. G. M., Buccheri, R., Caraveo, 

1319 P. A., Hermsen, W., Kanbach, G., Lebrun, F., Mayer-Hasselwander, H. A., & Strong, 

1320 A. W. 1985, Astron. & Astrophys., 146, 352 

1321 Pollock, A. M. T., Masnou, J. L., Bignami, G. F., Hermsen, W., Swanenburg, B. N., Kan- 

1322 bach, G., Lichti, G. G., & Wills, R. D. 1981, Astron. & Astrophys., 94, 116 

1323 Porciani, C, & Madau, P. 2001, Ap. J., 548, 522 

1324 Preece, R. D., Briggs, M. S., Giblin, T. W., Mallozzi, R. S., Pendleton, G. N., Paciesas, 

1325 W. S., & Band, D. L. 2002, Ap. J., 581, 1248 



-70- 



1326 Preece, R. D., Briggs, M. S., Mallozzi, R. S., Pendleton, G. N., Paciesas, W. S., & Band, 

1327 D. L. 1998, Ap. J. Lett., 506, L23 

1328 — . 2000, Ap. J. Supp., 126, 19 

1329 Rees, M. J., & Meszaros, P. 2005, Ap. J., 628, 847 

1330 Reyes, L. C. 2007, in American Institute of Physics Conference Series, Vol. 921, American 

1331 Institute of Physics Conference Series, ed. S. Ritz, P. Michelson, & C. A. Meegan, 

1332 359-360 

1333 Ryde, F. 2004, Ap. J., 614, 827 

1334 — . 2005, Ap. J. Lett., 625, L95 

1335 Scargle, J. D. 1998, Ap. J., 504, 405 

1336 Schaefer, B. E. 1999, Physical Review Letters, 82, 4964 

1337 Soderberg, A. M., Kulkarni, S. R., Nakar, E., Berger, E., Cameron, P. B., Fox, D. B., Frail, 

1338 D., Gal- Yam, A., Sari, R., Cenko, S. B., Kasliwal, M., Chevaher, R. A., Piran, T., 

1339 Price, P. A., Schmidt, B. P., Poolcy. G., Moon, D.-S., Pcnprasc, B. E., Ofck, E., Rau, 

1340 A., Gehrels, N., Nousek, J. A., Burrows, D. N., Persson, S. E., & McCarthy, P. J. 

1341 2006, Nature, 442, 1014 

1342 Sommer, M., Bertsch, D. L., Dingus, B. L., Fichtel, C. E., Fishman, G. J., Harding, A. K., 

1343 Hartman, R. C., Hunter, S. D., Hurley, K., Kanbach, G., Kniffen, D. A., Kouvehotou, 

1344 C, Lin, Y. C, Mattox, J. R., Maycr-Hassclwander, H. A., Michelson, P. F., von 

1345 Montigny, C, Nolan, P. L., Schneid, E., Sreekumar, P., & Thompson, D. J. 1994, Ap. 

1346 J. Lett., 422, L63 

1347 Stamatikos, M., & Band, D. L. 2006, in American Institute of Physics Conference Series, 

1348 Vol. 836, Gamma-Ray Bursts in the Swift Era, ed. S. S. Holt, N. Gehrels, & J. A. 

1349 Nousek, 599-604 

1350 Stamatikos, M., Sakamoto, T., & Band, D. L. 2008a, in American Institute of Physics 

1351 Conference Series, Vol. 1000, American Institute of Physics Conference Series, 133- 

1352 136 

1353 Stamatikos, M., Ukwatta, T. N., Sakamoto, T., Barthelmy, S. D., Norris, J. P., Gehrels, N., 

1354 & Dhuga, K. S. 2008b, in American Institute of Physics Conference Series, Vol. 1000, 

1355 American Institute of Physics Conference Series, 137-141 



1356 Stamatikos, M., et al. 2005, in International Cosmic Ray Conference, Vol. 4, International 

1357 Cosmic Ray Conference, 471 — h 

1358 Tajima, H., Bregeon, J., Chiang, J., & Thayer, G. 2008, GRB Coordinates Network, 8246, 1 

1359 Totani, T. 1998, Ap. J. Lett., 509, L81 

1360 Walker, K. C, Schaefer, B. E., & Fenimore, E. E. 2000, Ap. J., 537, 264 

1361 Wang, X. Y., Dai, Z. G., & Lu, T. 2001, Ap. J. Lett., 546, L33 

1362 Wang, X.-Y., Li, Z., & Meszaros, R 2006, Ap. J. Lett., 641, L89 

1363 Waxman, E. 1997, Ap. J. Lett., 485, L5+ 

1364 Waxman, E., & Bahcall, J. 1997, Physical Review Letters, 78, 2292 

1365 Wilks, S. S. 1938, Ann. Math. Stat., 9, 60 

1366 Xiao, L., & Schaefer, B. E. 2009, ArXiv e-prints 

1367 Yi, T., Liang, E., Qin, Y., & Lu, R. 2006, Mon. Not. RAS, 367, 1751 

1368 Zhang, B. 2007, Chinese Journal of Astronomy and Astrophysics, 7, 1 

1369 Zhang, B., & Meszaros, P. 2001, Ap. J., 559, 110 

1370 — . 2002, Ap. J., 581, 1236 



This preprint was prepared with the A AS IM^jX macros v5.2. 



